
%%% The Ambiguity of Fishing for Fun, October 2022
%%% This file computes and compares mean squared prediction errors for the four utility models considered

tic;
rng(4560); 

%Loading predicted data per state and year
MD2010=importdata('PredictedMD_2010.xls');MD2011=importdata('PredictedMD_2011.xls');MD2012=importdata('PredictedMD_2012.xls');MD2013=importdata('PredictedMD_2013.xls'); MD2014=importdata('PredictedMD_2014.xls'); MD2015=importdata('PredictedMD_2015.xls'); MD2016=importdata('PredictedMD_2016.xls'); MD2017=importdata('PredictedMD_2017.xls'); MD2018=importdata('PredictedMD_2018.xls'); MD2019=importdata('PredictedMD_2019.xls'); MD2020=importdata('PredictedMD_2020.xls');
MA2010=importdata('PredictedMA_2010.xls');MA2011=importdata('PredictedMA_2011.xls');MA2012=importdata('PredictedMA_2012.xls');MA2013=importdata('PredictedMA_2013.xls'); MA2014=importdata('PredictedMA_2014.xls'); MA2015=importdata('PredictedMA_2015.xls'); MA2016=importdata('PredictedMA_2016.xls'); MA2017=importdata('PredictedMA_2017.xls'); MA2018=importdata('PredictedMA_2018.xls'); MA2019=importdata('PredictedMA_2019.xls'); MA2020=importdata('PredictedMA_2020.xls');
NC2010=importdata('PredictedNC_2010.xls');NC2011=importdata('PredictedNC_2011.xls');NC2012=importdata('PredictedNC_2012.xls');NC2013=importdata('PredictedNC_2013.xls');NC2014=importdata('PredictedNC_2014.xls');NC2015=importdata('PredictedNC_2015.xls');NC2016=importdata('PredictedNC_2016.xls');NC2017=importdata('PredictedNC_2017.xls');NC2018=importdata('PredictedNC_2018.xls');NC2019=importdata('PredictedNC_2019.xls');NC2020=importdata('PredictedNC_2020.xls');
CT2010=importdata('PredictedCT_2010.xls');CT2011=importdata('PredictedCT_2011.xls');CT2012=importdata('PredictedCT_2012.xls');CT2013=importdata('PredictedCT_2013.xls');CT2014=importdata('PredictedCT_2014.xls');CT2015=importdata('PredictedCT_2015.xls');CT2016=importdata('PredictedCT_2016.xls');CT2017=importdata('PredictedCT_2017.xls'); CT2018=importdata('PredictedCT_2018.xls'); CT2019=importdata('PredictedCT_2019.xls'); CT2020=importdata('PredictedCT_2020.xls');
DE2010=importdata('PredictedDE_2010.xls');DE2011=importdata('PredictedDE_2011.xls');DE2012=importdata('PredictedDE_2012.xls');DE2013=importdata('PredictedDE_2013.xls');DE2014=importdata('PredictedDE_2014.xls');DE2015=importdata('PredictedDE_2015.xls');DE2016=importdata('PredictedDE_2016.xls');DE2017=importdata('PredictedDE_2017.xls');DE2018=importdata('PredictedDE_2018.xls');DE2019=importdata('PredictedDE_2019.xls');DE2020=importdata('PredictedDE_2020.xls');
VA2010=importdata('PredictedVA_2010.xls');VA2011=importdata('PredictedVA_2011.xls');VA2012=importdata('PredictedVA_2012.xls');VA2013=importdata('PredictedVA_2013.xls');VA2014=importdata('PredictedVA_2014.xls');VA2015=importdata('PredictedVA_2015.xls');VA2016=importdata('PredictedVA_2016.xls');VA2017=importdata('PredictedVA_2017.xls');VA2018=importdata('PredictedVA_2018.xls');VA2019=importdata('PredictedVA_2019.xls');VA2020=importdata('PredictedVA_2020.xls');
RI2010=importdata('PredictedRI_2010.xls');RI2011=importdata('PredictedRI_2011.xls');RI2012=importdata('PredictedRI_2012.xls');RI2013=importdata('PredictedRI_2013.xls');RI2014=importdata('PredictedRI_2014.xls');RI2015=importdata('PredictedRI_2015.xls');RI2016=importdata('PredictedRI_2016.xls');RI2017=importdata('PredictedRI_2017.xls');RI2018=importdata('PredictedRI_2018.xls');RI2019=importdata('PredictedRI_2019.xls');RI2020=importdata('PredictedRI_2020.xls');
NY2010=importdata('PredictedNY_2010.xls');NY2011=importdata('PredictedNY_2011.xls');NY2012=importdata('PredictedNY_2012.xls');NY2013=importdata('PredictedNY_2013.xls');NY2014=importdata('PredictedNY_2014.xls');NY2015=importdata('PredictedNY_2015.xls');NY2016=importdata('PredictedNY_2016.xls');NY2017=importdata('PredictedNY_2017.xls');NY2018=importdata('PredictedNY_2018.xls');NY2019=importdata('PredictedNY_2019.xls');NY2020=importdata('PredictedNY_2020.xls');
NJ2010=importdata('PredictedNJ_2010.xls');NJ2011=importdata('PredictedNJ_2011.xls');NJ2012=importdata('PredictedNJ_2012.xls');NJ2013=importdata('PredictedNJ_2013.xls');NJ2014=importdata('PredictedNJ_2014.xls');NJ2015=importdata('PredictedNJ_2015.xls');NJ2016=importdata('PredictedNJ_2016.xls');NJ2017=importdata('PredictedNJ_2017.xls');NJ2018=importdata('PredictedNJ_2018.xls');NJ2019=importdata('PredictedNJ_2019.xls');NJ2020=importdata('PredictedNJ_2020.xls');

%Loading actual data for each state from MRIP
%Trips
%Maryland
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'MARYLAND'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsMD10=NETrips.data(r10,2);
TripsMD11=NETrips.data(r11,2);
TripsMD12=NETrips.data(r12,2);
TripsMD13=NETrips.data(r13,2);
TripsMD14=NETrips.data(r14,2);
TripsMD15=NETrips.data(r15,2);
TripsMD16=NETrips.data(r16,2);
TripsMD17=NETrips.data(r17,2);
TripsMD18=NETrips.data(r18,2);
TripsMD19=NETrips.data(r19,2);
TripsMD20=NETrips.data(r20,2);

%Virginia
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'VIRGINIA'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsVA10=NETrips.data(r10,2);
TripsVA11=NETrips.data(r11,2);
TripsVA12=NETrips.data(r12,2);
TripsVA13=NETrips.data(r13,2);
TripsVA14=NETrips.data(r14,2);
TripsVA15=NETrips.data(r15,2);
TripsVA16=NETrips.data(r16,2);
TripsVA17=NETrips.data(r17,2);
TripsVA18=NETrips.data(r18,2);
TripsVA19=NETrips.data(r19,2);
TripsVA20=NETrips.data(r20,2);

%Massachusetts
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'MASSACHUSETTS'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsMA10=NETrips.data(r10,2);
TripsMA11=NETrips.data(r11,2);
TripsMA12=NETrips.data(r12,2);
TripsMA13=NETrips.data(r13,2);
TripsMA14=NETrips.data(r14,2);
TripsMA15=NETrips.data(r15,2);
TripsMA16=NETrips.data(r16,2);
TripsMA17=NETrips.data(r17,2);
TripsMA18=NETrips.data(r18,2);
TripsMA19=NETrips.data(r19,2);
TripsMA20=NETrips.data(r20,2);

%Delaware
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'DELAWARE'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,~]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsDE10=NETrips.data(r10,2);
TripsDE11=NETrips.data(r11,2);
TripsDE12=NETrips.data(r12,2);
TripsDE13=NETrips.data(r13,2);
TripsDE14=NETrips.data(r14,2);
TripsDE15=NETrips.data(r15,2);
TripsDE16=NETrips.data(r16,2);
TripsDE17=NETrips.data(r17,2);
TripsDE18=NETrips.data(r18,2);
TripsDE19=NETrips.data(r19,2);
TripsDE20=NETrips.data(r20,2);

%Connecticut
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'CONNECTICUT'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsCT10=NETrips.data(r10,2);
TripsCT11=NETrips.data(r11,2);
TripsCT12=NETrips.data(r12,2);
TripsCT13=NETrips.data(r13,2);
TripsCT14=NETrips.data(r14,2);
TripsCT15=NETrips.data(r15,2);
TripsCT16=NETrips.data(r16,2);
TripsCT17=NETrips.data(r17,2);
TripsCT18=NETrips.data(r18,2);
TripsCT19=NETrips.data(r19,2);
TripsCT20=NETrips.data(r20,2);

%North Carolina
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'NORTH CAROLINA'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsNC10=NETrips.data(r10,2);
TripsNC11=NETrips.data(r11,2);
TripsNC12=NETrips.data(r12,2);
TripsNC13=NETrips.data(r13,2);
TripsNC14=NETrips.data(r14,2);
TripsNC15=NETrips.data(r15,2);
TripsNC16=NETrips.data(r16,2);
TripsNC17=NETrips.data(r17,2);
TripsNC18=NETrips.data(r18,2);
TripsNC19=NETrips.data(r19,2);
TripsNC20=NETrips.data(r20,2);

%Rhode Island
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'RHODE ISLAND'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsRI10=NETrips.data(r10,2);
TripsRI11=NETrips.data(r11,2);
TripsRI12=NETrips.data(r12,2);
TripsRI13=NETrips.data(r13,2);
TripsRI14=NETrips.data(r14,2);
TripsRI15=NETrips.data(r15,2);
TripsRI16=NETrips.data(r16,2);
TripsRI17=NETrips.data(r17,2);
TripsRI18=NETrips.data(r18,2);
TripsRI19=NETrips.data(r19,2);
TripsRI20=NETrips.data(r20,2);

%New Jersey
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'NEW JERSEY'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsNJ10=NETrips.data(r10,2);
TripsNJ11=NETrips.data(r11,2);
TripsNJ12=NETrips.data(r12,2);
TripsNJ13=NETrips.data(r13,2);
TripsNJ14=NETrips.data(r14,2);
TripsNJ15=NETrips.data(r15,2);
TripsNJ16=NETrips.data(r16,2);
TripsNJ17=NETrips.data(r17,2);
TripsNJ18=NETrips.data(r18,2);
TripsNJ19=NETrips.data(r19,2);
TripsNJ20=NETrips.data(r20,2);

%New York
NETrips=importdata('mrip_directed_sf_trips_2009_2020.xlsx'); 
NETrips.textdata=NETrips.textdata(26:end,:);
[r9,c9]=find(contains(NETrips.textdata,'NEW YORK'));
NETrips.data=horzcat(NETrips.data(:,1),NETrips.data(:,3));
NETrips.data=NETrips.data(r9,:);
[r10,c10]=find(NETrips.data==2010);
[r11,c11]=find(NETrips.data==2011);
[r12,c12]=find(NETrips.data==2012);
[r13,c13]=find(NETrips.data==2013);
[r14,c14]=find(NETrips.data==2014);
[r15,c15]=find(NETrips.data==2015);
[r16,c16]=find(NETrips.data==2016);
[r17,c17]=find(NETrips.data==2017);
[r18,c18]=find(NETrips.data==2018);
[r19,c19]=find(NETrips.data==2019);
[r20,c20]=find(NETrips.data==2020);
TripsNY10=NETrips.data(r10,2);
TripsNY11=NETrips.data(r11,2);
TripsNY12=NETrips.data(r12,2);
TripsNY13=NETrips.data(r13,2);
TripsNY14=NETrips.data(r14,2);
TripsNY15=NETrips.data(r15,2);
TripsNY16=NETrips.data(r16,2);
TripsNY17=NETrips.data(r17,2);
TripsNY18=NETrips.data(r18,2);
TripsNY19=NETrips.data(r19,2);
TripsNY20=NETrips.data(r20,2);

ActualTrips10=(TripsNY10+TripsNJ10+TripsMD10+TripsNC10+TripsMA10+TripsVA10+TripsDE10+TripsRI10+TripsCT10);
ActualTrips11=(TripsNY11+TripsNJ11+TripsMD11+TripsNC11+TripsMA11+TripsVA11+TripsDE11+TripsRI11+TripsCT11);
ActualTrips12=(TripsNY12+TripsNJ12+TripsMD12+TripsNC12+TripsMA12+TripsVA12+TripsDE12+TripsRI12+TripsCT12);
ActualTrips13=(TripsNY13+TripsNJ13+TripsMD13+TripsNC13+TripsMA13+TripsVA13+TripsDE13+TripsRI13+TripsCT13);
ActualTrips14=(TripsNY14+TripsNJ14+TripsMD14+TripsNC14+TripsMA14+TripsVA14+TripsDE14+TripsRI14+TripsCT14);
ActualTrips15=(TripsNY15+TripsNJ15+TripsMD15+TripsNC15+TripsMA15+TripsVA15+TripsDE15+TripsRI15+TripsCT15);
ActualTrips16=(TripsNY16+TripsNJ16+TripsMD16+TripsNC16+TripsMA16+TripsVA16+TripsDE16+TripsRI16+TripsCT16);
ActualTrips17=(TripsNY17+TripsNJ17+TripsMD17+TripsNC17+TripsMA17+TripsVA17+TripsDE17+TripsRI17+TripsCT17);
ActualTrips18=(TripsNY18+TripsNJ18+TripsMD18+TripsNC18+TripsMA18+TripsVA18+TripsDE18+TripsRI18+TripsCT18);
ActualTrips19=(TripsNY19+TripsNJ19+TripsMD19+TripsNC19+TripsMA19+TripsVA19+TripsDE19+TripsRI19+TripsCT19);
ActualTrips20=(TripsNY20+TripsNJ20+TripsMD20+TripsNC20+TripsMA20+TripsVA20+TripsDE20+TripsRI20+TripsCT20);

Mean_ActualTrips=(ActualTrips10+ActualTrips11+ActualTrips12+ActualTrips13+ActualTrips14+ActualTrips15+ActualTrips16+ActualTrips17+ActualTrips18+ActualTrips19+ActualTrips20)/1000/11;

%Harvest (numbers of fish)
%New York
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'NEW YORK'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishNY10=NEFish.data(r10,2);
FishNY11=NEFish.data(r11,2);
FishNY12=NEFish.data(r12,2);
FishNY13=NEFish.data(r13,2);
FishNY14=NEFish.data(r14,2);
FishNY15=NEFish.data(r15,2);
FishNY16=NEFish.data(r16,2);
FishNY17=NEFish.data(r17,2);
FishNY18=NEFish.data(r18,2);
FishNY19=NEFish.data(r19,2);
FishNY20=NEFish.data(r20,2);

%New Jersey
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'NEW JERSEY'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishNJ10=NEFish.data(r10,2);
FishNJ11=NEFish.data(r11,2);
FishNJ12=NEFish.data(r12,2);
FishNJ13=NEFish.data(r13,2);
FishNJ14=NEFish.data(r14,2);
FishNJ15=NEFish.data(r15,2);
FishNJ16=NEFish.data(r16,2);
FishNJ17=NEFish.data(r17,2);
FishNJ18=NEFish.data(r18,2);
FishNJ19=NEFish.data(r19,2);
FishNJ20=NEFish.data(r20,2);

%Maryland
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'MARYLAND'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishMD10=NEFish.data(r10,2);
FishMD11=NEFish.data(r11,2);
FishMD12=NEFish.data(r12,2);
FishMD13=NEFish.data(r13,2);
FishMD14=NEFish.data(r14,2);
FishMD15=NEFish.data(r15,2);
FishMD16=NEFish.data(r16,2);
FishMD17=NEFish.data(r17,2);
FishMD18=NEFish.data(r18,2);
FishMD19=NEFish.data(r19,2);
FishMD20=NEFish.data(r20,2);

%Delaware
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'DELAWARE'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishDE10=NEFish.data(r10,2);
FishDE11=NEFish.data(r11,2);
FishDE12=NEFish.data(r12,2);
FishDE13=NEFish.data(r13,2);
FishDE14=NEFish.data(r14,2);
FishDE15=NEFish.data(r15,2);
FishDE16=NEFish.data(r16,2);
FishDE17=NEFish.data(r17,2);
FishDE18=NEFish.data(r18,2);
FishDE19=NEFish.data(r19,2);
FishDE20=NEFish.data(r20,2);

%North Carolina
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'NORTH CAROLINA'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishNC10=NEFish.data(r10,2);
FishNC11=NEFish.data(r11,2);
FishNC12=NEFish.data(r12,2);
FishNC13=NEFish.data(r13,2);
FishNC14=NEFish.data(r14,2);
FishNC15=NEFish.data(r15,2);
FishNC16=NEFish.data(r16,2);
FishNC17=NEFish.data(r17,2);
FishNC18=NEFish.data(r18,2);
FishNC19=NEFish.data(r19,2);
FishNC20=NEFish.data(r20,2);

%Virginia
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'VIRGINIA'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishVA10=NEFish.data(r10,2);
FishVA11=NEFish.data(r11,2);
FishVA12=NEFish.data(r12,2);
FishVA13=NEFish.data(r13,2);
FishVA14=NEFish.data(r14,2);
FishVA15=NEFish.data(r15,2);
FishVA16=NEFish.data(r16,2);
FishVA17=NEFish.data(r17,2);
FishVA18=NEFish.data(r18,2);
FishVA19=NEFish.data(r19,2);
FishVA20=NEFish.data(r20,2);

%Massachusetts
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'MASSACHUSETTS'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishMA10=NEFish.data(r10,2);
FishMA11=NEFish.data(r11,2);
FishMA12=NEFish.data(r12,2);
FishMA13=NEFish.data(r13,2);
FishMA14=NEFish.data(r14,2);
FishMA15=NEFish.data(r15,2);
FishMA16=NEFish.data(r16,2);
FishMA17=NEFish.data(r17,2);
FishMA18=NEFish.data(r18,2);
FishMA19=NEFish.data(r19,2);
FishMA20=NEFish.data(r20,2);

%Connecticut
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'CONNECTICUT'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishCT10=NEFish.data(r10,2);
FishCT11=NEFish.data(r11,2);
FishCT12=NEFish.data(r12,2);
FishCT13=NEFish.data(r13,2);
FishCT14=NEFish.data(r14,2);
FishCT15=NEFish.data(r15,2);
FishCT16=NEFish.data(r16,2);
FishCT17=NEFish.data(r17,2);
FishCT18=NEFish.data(r18,2);
FishCT19=NEFish.data(r19,2);
FishCT20=NEFish.data(r20,2);

%Rhode Island
NEFish=importdata('mrip_SUMMER_FLOUNDER_catch_series.xls'); 
NEFish.textdata=NEFish.textdata(30:end,:);
[r9,c9]=find(contains(NEFish.textdata,'RHODE ISLAND'));
NEFish.data=horzcat(NEFish.data(:,1),NEFish.data(:,4));
NEFish.data=NEFish.data(r9,:);
[r10,c10]=find(NEFish.data==2010);
[r11,c11]=find(NEFish.data==2011);
[r12,c12]=find(NEFish.data==2012);
[r13,c13]=find(NEFish.data==2013);
[r14,c14]=find(NEFish.data==2014);
[r15,c15]=find(NEFish.data==2015);
[r16,c16]=find(NEFish.data==2016);
[r17,c17]=find(NEFish.data==2017);
[r18,c18]=find(NEFish.data==2018);
[r19,c19]=find(NEFish.data==2019);
[r20,c20]=find(NEFish.data==2020);
FishRI10=NEFish.data(r10,2);
FishRI11=NEFish.data(r11,2);
FishRI12=NEFish.data(r12,2);
FishRI13=NEFish.data(r13,2);
FishRI14=NEFish.data(r14,2);
FishRI15=NEFish.data(r15,2);
FishRI16=NEFish.data(r16,2);
FishRI17=NEFish.data(r17,2);
FishRI18=NEFish.data(r18,2);
FishRI19=NEFish.data(r19,2);
FishRI20=NEFish.data(r20,2);

ActualFish10=(FishNY10+FishNJ10+FishMD10+FishNC10+FishMA10+FishVA10+FishDE10+FishRI10+FishCT10);
ActualFish11=(FishNY11+FishNJ11+FishMD11+FishNC11+FishMA11+FishVA11+FishDE11+FishRI11+FishCT11);
ActualFish12=(FishNY12+FishNJ12+FishMD12+FishNC12+FishMA12+FishVA12+FishDE12+FishRI12+FishCT12);
ActualFish13=(FishNY13+FishNJ13+FishMD13+FishNC13+FishMA13+FishVA13+FishDE13+FishRI13+FishCT13);
ActualFish14=(FishNY14+FishNJ14+FishMD14+FishNC14+FishMA14+FishVA14+FishDE14+FishRI14+FishCT14);
ActualFish15=(FishNY15+FishNJ15+FishMD15+FishNC15+FishMA15+FishVA15+FishDE15+FishRI15+FishCT15);
ActualFish16=(FishNY16+FishNJ16+FishMD16+FishNC16+FishMA16+FishVA16+FishDE16+FishRI16+FishCT16);
ActualFish17=(FishNY17+FishNJ17+FishMD17+FishNC17+FishMA17+FishVA17+FishDE17+FishRI17+FishCT17);
ActualFish18=(FishNY18+FishNJ18+FishMD18+FishNC18+FishMA18+FishVA18+FishDE18+FishRI18+FishCT18);
ActualFish19=(FishNY19+FishNJ19+FishMD19+FishNC19+FishMA19+FishVA19+FishDE19+FishRI19+FishCT19);
ActualFish20=(FishNY20+FishNJ20+FishMD20+FishNC20+FishMA20+FishVA20+FishDE20+FishRI20+FishCT20);

Mean_ActualFish=(ActualFish10+ActualFish11+ActualFish12+ActualFish13+ActualFish14+ActualFish15+ActualFish16+ActualFish17+ActualFish18+ActualFish19+ActualFish20)/1000/11;


%Predicted data
%Trips
PredictedTrips10A=MD2010.data.Sheet1(:,1)+MA2010.data.Sheet1(:,1)+DE2010.data.Sheet1(:,1)+NC2010.data.Sheet1(:,1)+CT2010.data.Sheet1(:,1)+RI2010.data.Sheet1(:,1)+VA2010.data.Sheet1(:,1)+NJ2010.data.Sheet1(:,1)+NY2010.data.Sheet1(:,1);
PredictedTrips10M=MD2010.data.Sheet2(:,1)+MA2010.data.Sheet2(:,1)+DE2010.data.Sheet2(:,1)+NC2010.data.Sheet2(:,1)+CT2010.data.Sheet2(:,1)+RI2010.data.Sheet2(:,1)+VA2010.data.Sheet2(:,1)+NJ2010.data.Sheet2(:,1)+NY2010.data.Sheet2(:,1);
PredictedTrips10L=MD2010.data.Sheet3(:,1)+MA2010.data.Sheet3(:,1)+DE2010.data.Sheet3(:,1)+NC2010.data.Sheet3(:,1)+CT2010.data.Sheet3(:,1)+RI2010.data.Sheet3(:,1)+VA2010.data.Sheet3(:,1)+NJ2010.data.Sheet3(:,1)+NY2010.data.Sheet3(:,1);
PredictedTrips10C=MD2010.data.Sheet4(:,1)+MA2010.data.Sheet4(:,1)+DE2010.data.Sheet4(:,1)+NC2010.data.Sheet4(:,1)+CT2010.data.Sheet4(:,1)+RI2010.data.Sheet4(:,1)+VA2010.data.Sheet4(:,1)+NJ2010.data.Sheet4(:,1)+NY2010.data.Sheet4(:,1);

PredictedTrips11A=MD2011.data.Sheet1(:,1)+MA2011.data.Sheet1(:,1)+DE2011.data.Sheet1(:,1)+NC2011.data.Sheet1(:,1)+CT2011.data.Sheet1(:,1)+RI2011.data.Sheet1(:,1)+VA2011.data.Sheet1(:,1)+NJ2011.data.Sheet1(:,1)+NY2011.data.Sheet1(:,1);
PredictedTrips11M=MD2011.data.Sheet2(:,1)+MA2011.data.Sheet2(:,1)+DE2011.data.Sheet2(:,1)+NC2011.data.Sheet2(:,1)+CT2011.data.Sheet2(:,1)+RI2011.data.Sheet2(:,1)+VA2011.data.Sheet2(:,1)+NJ2011.data.Sheet2(:,1)+NY2011.data.Sheet2(:,1);
PredictedTrips11L=MD2011.data.Sheet3(:,1)+MA2011.data.Sheet3(:,1)+DE2011.data.Sheet3(:,1)+NC2011.data.Sheet3(:,1)+CT2011.data.Sheet3(:,1)+RI2011.data.Sheet3(:,1)+VA2011.data.Sheet3(:,1)+NJ2011.data.Sheet3(:,1)+NY2011.data.Sheet3(:,1);
PredictedTrips11C=MD2011.data.Sheet4(:,1)+MA2011.data.Sheet4(:,1)+DE2011.data.Sheet4(:,1)+NC2011.data.Sheet4(:,1)+CT2011.data.Sheet4(:,1)+RI2011.data.Sheet4(:,1)+VA2011.data.Sheet4(:,1)+NJ2011.data.Sheet4(:,1)+NY2011.data.Sheet4(:,1);

PredictedTrips12A=MD2012.data.Sheet1(:,1)+MA2012.data.Sheet1(:,1)+DE2012.data.Sheet1(:,1)+NC2012.data.Sheet1(:,1)+CT2012.data.Sheet1(:,1)+RI2012.data.Sheet1(:,1)+VA2012.data.Sheet1(:,1)+NJ2012.data.Sheet1(:,1)+NY2012.data.Sheet1(:,1);
PredictedTrips12M=MD2012.data.Sheet2(:,1)+MA2012.data.Sheet2(:,1)+DE2012.data.Sheet2(:,1)+NC2012.data.Sheet2(:,1)+CT2012.data.Sheet2(:,1)+RI2012.data.Sheet2(:,1)+VA2012.data.Sheet2(:,1)+NJ2012.data.Sheet2(:,1)+NY2012.data.Sheet2(:,1);
PredictedTrips12L=MD2012.data.Sheet3(:,1)+MA2012.data.Sheet3(:,1)+DE2012.data.Sheet3(:,1)+NC2012.data.Sheet3(:,1)+CT2012.data.Sheet3(:,1)+RI2012.data.Sheet3(:,1)+VA2012.data.Sheet3(:,1)+NJ2012.data.Sheet3(:,1)+NY2012.data.Sheet3(:,1);
PredictedTrips12C=MD2012.data.Sheet4(:,1)+MA2012.data.Sheet4(:,1)+DE2012.data.Sheet4(:,1)+NC2012.data.Sheet4(:,1)+CT2012.data.Sheet4(:,1)+RI2012.data.Sheet4(:,1)+VA2012.data.Sheet4(:,1)+NJ2012.data.Sheet4(:,1)+NY2012.data.Sheet4(:,1);

PredictedTrips13A=MD2013.data.Sheet1(:,1)+MA2013.data.Sheet1(:,1)+DE2013.data.Sheet1(:,1)+NC2013.data.Sheet1(:,1)+CT2013.data.Sheet1(:,1)+RI2013.data.Sheet1(:,1)+VA2013.data.Sheet1(:,1)+NJ2013.data.Sheet1(:,1)+NY2013.data.Sheet1(:,1);
PredictedTrips13M=MD2013.data.Sheet2(:,1)+MA2013.data.Sheet2(:,1)+DE2013.data.Sheet2(:,1)+NC2013.data.Sheet2(:,1)+CT2013.data.Sheet2(:,1)+RI2013.data.Sheet2(:,1)+VA2013.data.Sheet2(:,1)+NJ2013.data.Sheet2(:,1)+NY2013.data.Sheet2(:,1);
PredictedTrips13L=MD2013.data.Sheet3(:,1)+MA2013.data.Sheet3(:,1)+DE2013.data.Sheet3(:,1)+NC2013.data.Sheet3(:,1)+CT2013.data.Sheet3(:,1)+RI2013.data.Sheet3(:,1)+VA2013.data.Sheet3(:,1)+NJ2013.data.Sheet3(:,1)+NY2013.data.Sheet3(:,1);
PredictedTrips13C=MD2013.data.Sheet4(:,1)+MA2013.data.Sheet4(:,1)+DE2013.data.Sheet4(:,1)+NC2013.data.Sheet4(:,1)+CT2013.data.Sheet4(:,1)+RI2013.data.Sheet4(:,1)+VA2013.data.Sheet4(:,1)+NJ2013.data.Sheet4(:,1)+NY2013.data.Sheet4(:,1);

PredictedTrips14A=MD2014.data.Sheet1(:,1)+MA2014.data.Sheet1(:,1)+DE2014.data.Sheet1(:,1)+NC2014.data.Sheet1(:,1)+CT2014.data.Sheet1(:,1)+RI2014.data.Sheet1(:,1)+VA2014.data.Sheet1(:,1)+NJ2014.data.Sheet1(:,1)+NY2014.data.Sheet1(:,1);
PredictedTrips14M=MD2014.data.Sheet2(:,1)+MA2014.data.Sheet2(:,1)+DE2014.data.Sheet2(:,1)+NC2014.data.Sheet2(:,1)+CT2014.data.Sheet2(:,1)+RI2014.data.Sheet2(:,1)+VA2014.data.Sheet2(:,1)+NJ2014.data.Sheet2(:,1)+NY2014.data.Sheet2(:,1);
PredictedTrips14L=MD2014.data.Sheet3(:,1)+MA2014.data.Sheet3(:,1)+DE2014.data.Sheet3(:,1)+NC2014.data.Sheet3(:,1)+CT2014.data.Sheet3(:,1)+RI2014.data.Sheet3(:,1)+VA2014.data.Sheet3(:,1)+NJ2014.data.Sheet3(:,1)+NY2014.data.Sheet3(:,1);
PredictedTrips14C=MD2014.data.Sheet4(:,1)+MA2014.data.Sheet4(:,1)+DE2014.data.Sheet4(:,1)+NC2014.data.Sheet4(:,1)+CT2014.data.Sheet4(:,1)+RI2014.data.Sheet4(:,1)+VA2014.data.Sheet4(:,1)+NJ2014.data.Sheet4(:,1)+NY2014.data.Sheet4(:,1);

PredictedTrips15A=MD2015.data.Sheet1(:,1)+MA2015.data.Sheet1(:,1)+DE2015.data.Sheet1(:,1)+NC2015.data.Sheet1(:,1)+CT2015.data.Sheet1(:,1)+RI2015.data.Sheet1(:,1)+VA2015.data.Sheet1(:,1)+NJ2015.data.Sheet1(:,1)+NY2015.data.Sheet1(:,1);
PredictedTrips15M=MD2015.data.Sheet2(:,1)+MA2015.data.Sheet2(:,1)+DE2015.data.Sheet2(:,1)+NC2015.data.Sheet2(:,1)+CT2015.data.Sheet2(:,1)+RI2015.data.Sheet2(:,1)+VA2015.data.Sheet2(:,1)+NJ2015.data.Sheet2(:,1)+NY2015.data.Sheet2(:,1);
PredictedTrips15L=MD2015.data.Sheet3(:,1)+MA2015.data.Sheet3(:,1)+DE2015.data.Sheet3(:,1)+NC2015.data.Sheet3(:,1)+CT2015.data.Sheet3(:,1)+RI2015.data.Sheet3(:,1)+VA2015.data.Sheet3(:,1)+NJ2015.data.Sheet3(:,1)+NY2015.data.Sheet3(:,1);
PredictedTrips15C=MD2015.data.Sheet4(:,1)+MA2015.data.Sheet4(:,1)+DE2015.data.Sheet4(:,1)+NC2015.data.Sheet4(:,1)+CT2015.data.Sheet4(:,1)+RI2015.data.Sheet4(:,1)+VA2015.data.Sheet4(:,1)+NJ2015.data.Sheet4(:,1)+NY2015.data.Sheet4(:,1);

PredictedTrips16A=MD2016.data.Sheet1(:,1)+MA2016.data.Sheet1(:,1)+DE2016.data.Sheet1(:,1)+NC2016.data.Sheet1(:,1)+CT2016.data.Sheet1(:,1)+RI2016.data.Sheet1(:,1)+VA2016.data.Sheet1(:,1)+NJ2016.data.Sheet1(:,1)+NY2016.data.Sheet1(:,1);
PredictedTrips16M=MD2016.data.Sheet2(:,1)+MA2016.data.Sheet2(:,1)+DE2016.data.Sheet2(:,1)+NC2016.data.Sheet2(:,1)+CT2016.data.Sheet2(:,1)+RI2016.data.Sheet2(:,1)+VA2016.data.Sheet2(:,1)+NJ2016.data.Sheet2(:,1)+NY2016.data.Sheet2(:,1);
PredictedTrips16L=MD2016.data.Sheet3(:,1)+MA2016.data.Sheet3(:,1)+DE2016.data.Sheet3(:,1)+NC2016.data.Sheet3(:,1)+CT2016.data.Sheet3(:,1)+RI2016.data.Sheet3(:,1)+VA2016.data.Sheet3(:,1)+NJ2016.data.Sheet3(:,1)+NY2016.data.Sheet3(:,1);
PredictedTrips16C=MD2016.data.Sheet4(:,1)+MA2016.data.Sheet4(:,1)+DE2016.data.Sheet4(:,1)+NC2016.data.Sheet4(:,1)+CT2016.data.Sheet4(:,1)+RI2016.data.Sheet4(:,1)+VA2016.data.Sheet4(:,1)+NJ2016.data.Sheet4(:,1)+NY2016.data.Sheet4(:,1);

PredictedTrips17A=MD2017.data.Sheet1(:,1)+MA2017.data.Sheet1(:,1)+DE2017.data.Sheet1(:,1)+NC2017.data.Sheet1(:,1)+CT2017.data.Sheet1(:,1)+RI2017.data.Sheet1(:,1)+VA2017.data.Sheet1(:,1)+NJ2017.data.Sheet1(:,1)+NY2017.data.Sheet1(:,1);
PredictedTrips17M=MD2017.data.Sheet2(:,1)+MA2017.data.Sheet2(:,1)+DE2017.data.Sheet2(:,1)+NC2017.data.Sheet2(:,1)+CT2017.data.Sheet2(:,1)+RI2017.data.Sheet2(:,1)+VA2017.data.Sheet2(:,1)+NJ2017.data.Sheet2(:,1)+NY2017.data.Sheet2(:,1);
PredictedTrips17L=MD2017.data.Sheet3(:,1)+MA2017.data.Sheet3(:,1)+DE2017.data.Sheet3(:,1)+NC2017.data.Sheet3(:,1)+CT2017.data.Sheet3(:,1)+RI2017.data.Sheet3(:,1)+VA2017.data.Sheet3(:,1)+NJ2017.data.Sheet3(:,1)+NY2017.data.Sheet3(:,1);
PredictedTrips17C=MD2017.data.Sheet4(:,1)+MA2017.data.Sheet4(:,1)+DE2017.data.Sheet4(:,1)+NC2017.data.Sheet4(:,1)+CT2017.data.Sheet4(:,1)+RI2017.data.Sheet4(:,1)+VA2017.data.Sheet4(:,1)+NJ2017.data.Sheet4(:,1)+NY2017.data.Sheet4(:,1);

PredictedTrips18A=MD2018.data.Sheet1(:,1)+MA2018.data.Sheet1(:,1)+DE2018.data.Sheet1(:,1)+NC2018.data.Sheet1(:,1)+CT2018.data.Sheet1(:,1)+RI2018.data.Sheet1(:,1)+VA2018.data.Sheet1(:,1)+NJ2018.data.Sheet1(:,1)+NY2018.data.Sheet1(:,1);
PredictedTrips18M=MD2018.data.Sheet2(:,1)+MA2018.data.Sheet2(:,1)+DE2018.data.Sheet2(:,1)+NC2018.data.Sheet2(:,1)+CT2018.data.Sheet2(:,1)+RI2018.data.Sheet2(:,1)+VA2018.data.Sheet2(:,1)+NJ2018.data.Sheet2(:,1)+NY2018.data.Sheet2(:,1);
PredictedTrips18L=MD2018.data.Sheet3(:,1)+MA2018.data.Sheet3(:,1)+DE2018.data.Sheet3(:,1)+NC2018.data.Sheet3(:,1)+CT2018.data.Sheet3(:,1)+RI2018.data.Sheet3(:,1)+VA2018.data.Sheet3(:,1)+NJ2018.data.Sheet3(:,1)+NY2018.data.Sheet3(:,1);
PredictedTrips18C=MD2018.data.Sheet4(:,1)+MA2018.data.Sheet4(:,1)+DE2018.data.Sheet4(:,1)+NC2018.data.Sheet4(:,1)+CT2018.data.Sheet4(:,1)+RI2018.data.Sheet4(:,1)+VA2018.data.Sheet4(:,1)+NJ2018.data.Sheet4(:,1)+NY2018.data.Sheet4(:,1);

PredictedTrips19A=MD2019.data.Sheet1(:,1)+MA2019.data.Sheet1(:,1)+DE2019.data.Sheet1(:,1)+NC2019.data.Sheet1(:,1)+CT2019.data.Sheet1(:,1)+RI2019.data.Sheet1(:,1)+VA2019.data.Sheet1(:,1)+NJ2019.data.Sheet1(:,1)+NY2019.data.Sheet1(:,1);
PredictedTrips19M=MD2019.data.Sheet2(:,1)+MA2019.data.Sheet2(:,1)+DE2019.data.Sheet2(:,1)+NC2019.data.Sheet2(:,1)+CT2019.data.Sheet2(:,1)+RI2019.data.Sheet2(:,1)+VA2019.data.Sheet2(:,1)+NJ2019.data.Sheet2(:,1)+NY2019.data.Sheet2(:,1);
PredictedTrips19L=MD2019.data.Sheet3(:,1)+MA2019.data.Sheet3(:,1)+DE2019.data.Sheet3(:,1)+NC2019.data.Sheet3(:,1)+CT2019.data.Sheet3(:,1)+RI2019.data.Sheet3(:,1)+VA2019.data.Sheet3(:,1)+NJ2019.data.Sheet3(:,1)+NY2019.data.Sheet3(:,1);
PredictedTrips19C=MD2019.data.Sheet4(:,1)+MA2019.data.Sheet4(:,1)+DE2019.data.Sheet4(:,1)+NC2019.data.Sheet4(:,1)+CT2019.data.Sheet4(:,1)+RI2019.data.Sheet4(:,1)+VA2019.data.Sheet4(:,1)+NJ2019.data.Sheet4(:,1)+NY2019.data.Sheet4(:,1);

PredictedTrips20A=MD2020.data.Sheet1(:,1)+MA2020.data.Sheet1(:,1)+DE2020.data.Sheet1(:,1)+NC2020.data.Sheet1(:,1)+CT2020.data.Sheet1(:,1)+RI2020.data.Sheet1(:,1)+VA2020.data.Sheet1(:,1)+NJ2020.data.Sheet1(:,1)+NY2020.data.Sheet1(:,1);
PredictedTrips20M=MD2020.data.Sheet2(:,1)+MA2020.data.Sheet2(:,1)+DE2020.data.Sheet2(:,1)+NC2020.data.Sheet2(:,1)+CT2020.data.Sheet2(:,1)+RI2020.data.Sheet2(:,1)+VA2020.data.Sheet2(:,1)+NJ2020.data.Sheet2(:,1)+NY2020.data.Sheet2(:,1);
PredictedTrips20L=MD2020.data.Sheet3(:,1)+MA2020.data.Sheet3(:,1)+DE2020.data.Sheet3(:,1)+NC2020.data.Sheet3(:,1)+CT2020.data.Sheet3(:,1)+RI2020.data.Sheet3(:,1)+VA2020.data.Sheet3(:,1)+NJ2020.data.Sheet3(:,1)+NY2020.data.Sheet3(:,1);
PredictedTrips20C=MD2020.data.Sheet4(:,1)+MA2020.data.Sheet4(:,1)+DE2020.data.Sheet4(:,1)+NC2020.data.Sheet4(:,1)+CT2020.data.Sheet4(:,1)+RI2020.data.Sheet4(:,1)+VA2020.data.Sheet4(:,1)+NJ2020.data.Sheet4(:,1)+NY2020.data.Sheet4(:,1);

%Landings:
PredictedFish10A=MD2010.data.Sheet1(:,4)+MA2010.data.Sheet1(:,4)+DE2010.data.Sheet1(:,4)+NC2010.data.Sheet1(:,4)+CT2010.data.Sheet1(:,4)+RI2010.data.Sheet1(:,4)+VA2010.data.Sheet1(:,4)+NJ2010.data.Sheet1(:,4)+NY2010.data.Sheet1(:,4);
PredictedFish10M=MD2010.data.Sheet2(:,4)+MA2010.data.Sheet2(:,4)+DE2010.data.Sheet2(:,4)+NC2010.data.Sheet2(:,4)+CT2010.data.Sheet2(:,4)+RI2010.data.Sheet2(:,4)+VA2010.data.Sheet2(:,4)+NJ2010.data.Sheet2(:,4)+NY2010.data.Sheet2(:,4);
PredictedFish10L=MD2010.data.Sheet3(:,4)+MA2010.data.Sheet3(:,4)+DE2010.data.Sheet3(:,4)+NC2010.data.Sheet3(:,4)+CT2010.data.Sheet3(:,4)+RI2010.data.Sheet3(:,4)+VA2010.data.Sheet3(:,4)+NJ2010.data.Sheet3(:,4)+NY2010.data.Sheet3(:,4);
PredictedFish10C=MD2010.data.Sheet4(:,4)+MA2010.data.Sheet4(:,4)+DE2010.data.Sheet4(:,4)+NC2010.data.Sheet4(:,4)+CT2010.data.Sheet4(:,4)+RI2010.data.Sheet4(:,4)+VA2010.data.Sheet4(:,4)+NJ2010.data.Sheet4(:,4)+NY2010.data.Sheet4(:,4);

PredictedFish11A=MD2011.data.Sheet1(:,4)+MA2011.data.Sheet1(:,4)+DE2011.data.Sheet1(:,4)+NC2011.data.Sheet1(:,4)+CT2011.data.Sheet1(:,4)+RI2011.data.Sheet1(:,4)+VA2011.data.Sheet1(:,4)+NJ2011.data.Sheet1(:,4)+NY2011.data.Sheet1(:,4);
PredictedFish11M=MD2011.data.Sheet2(:,4)+MA2011.data.Sheet2(:,4)+DE2011.data.Sheet2(:,4)+NC2011.data.Sheet2(:,4)+CT2011.data.Sheet2(:,4)+RI2011.data.Sheet2(:,4)+VA2011.data.Sheet2(:,4)+NJ2011.data.Sheet2(:,4)+NY2011.data.Sheet2(:,4);
PredictedFish11L=MD2011.data.Sheet3(:,4)+MA2011.data.Sheet3(:,4)+DE2011.data.Sheet3(:,4)+NC2011.data.Sheet3(:,4)+CT2011.data.Sheet3(:,4)+RI2011.data.Sheet3(:,4)+VA2011.data.Sheet3(:,4)+NJ2011.data.Sheet3(:,4)+NY2011.data.Sheet3(:,4);
PredictedFish11C=MD2011.data.Sheet4(:,4)+MA2011.data.Sheet4(:,4)+DE2011.data.Sheet4(:,4)+NC2011.data.Sheet4(:,4)+CT2011.data.Sheet4(:,4)+RI2011.data.Sheet4(:,4)+VA2011.data.Sheet4(:,4)+NJ2011.data.Sheet4(:,4)+NY2011.data.Sheet4(:,4);

PredictedFish12A=MD2012.data.Sheet1(:,4)+MA2012.data.Sheet1(:,4)+DE2012.data.Sheet1(:,4)+NC2012.data.Sheet1(:,4)+CT2012.data.Sheet1(:,4)+RI2012.data.Sheet1(:,4)+VA2012.data.Sheet1(:,4)+NJ2012.data.Sheet1(:,4)+NY2012.data.Sheet1(:,4);
PredictedFish12M=MD2012.data.Sheet2(:,4)+MA2012.data.Sheet2(:,4)+DE2012.data.Sheet2(:,4)+NC2012.data.Sheet2(:,4)+CT2012.data.Sheet2(:,4)+RI2012.data.Sheet2(:,4)+VA2012.data.Sheet2(:,4)+NJ2012.data.Sheet2(:,4)+NY2012.data.Sheet2(:,4);
PredictedFish12L=MD2012.data.Sheet3(:,4)+MA2012.data.Sheet3(:,4)+DE2012.data.Sheet3(:,4)+NC2012.data.Sheet3(:,4)+CT2012.data.Sheet3(:,4)+RI2012.data.Sheet3(:,4)+VA2012.data.Sheet3(:,4)+NJ2012.data.Sheet3(:,4)+NY2012.data.Sheet3(:,4);
PredictedFish12C=MD2012.data.Sheet4(:,4)+MA2012.data.Sheet4(:,4)+DE2012.data.Sheet4(:,4)+NC2012.data.Sheet4(:,4)+CT2012.data.Sheet4(:,4)+RI2012.data.Sheet4(:,4)+VA2012.data.Sheet4(:,4)+NJ2012.data.Sheet4(:,4)+NY2012.data.Sheet4(:,4);

PredictedFish13A=MD2013.data.Sheet1(:,4)+MA2013.data.Sheet1(:,4)+DE2013.data.Sheet1(:,4)+NC2013.data.Sheet1(:,4)+CT2013.data.Sheet1(:,4)+RI2013.data.Sheet1(:,4)+VA2013.data.Sheet1(:,4)+NJ2013.data.Sheet1(:,4)+NY2013.data.Sheet1(:,4);
PredictedFish13M=MD2013.data.Sheet2(:,4)+MA2013.data.Sheet2(:,4)+DE2013.data.Sheet2(:,4)+NC2013.data.Sheet2(:,4)+CT2013.data.Sheet2(:,4)+RI2013.data.Sheet2(:,4)+VA2013.data.Sheet2(:,4)+NJ2013.data.Sheet2(:,4)+NY2013.data.Sheet2(:,4);
PredictedFish13L=MD2013.data.Sheet3(:,4)+MA2013.data.Sheet3(:,4)+DE2013.data.Sheet3(:,4)+NC2013.data.Sheet3(:,4)+CT2013.data.Sheet3(:,4)+RI2013.data.Sheet3(:,4)+VA2013.data.Sheet3(:,4)+NJ2013.data.Sheet3(:,4)+NY2013.data.Sheet3(:,4);
PredictedFish13C=MD2013.data.Sheet4(:,4)+MA2013.data.Sheet4(:,4)+DE2013.data.Sheet4(:,4)+NC2013.data.Sheet4(:,4)+CT2013.data.Sheet4(:,4)+RI2013.data.Sheet4(:,4)+VA2013.data.Sheet4(:,4)+NJ2013.data.Sheet4(:,4)+NY2013.data.Sheet4(:,4);

PredictedFish14A=MD2014.data.Sheet1(:,4)+MA2014.data.Sheet1(:,4)+DE2014.data.Sheet1(:,4)+NC2014.data.Sheet1(:,4)+CT2014.data.Sheet1(:,4)+RI2014.data.Sheet1(:,4)+VA2014.data.Sheet1(:,4)+NJ2014.data.Sheet1(:,4)+NY2014.data.Sheet1(:,4);
PredictedFish14M=MD2014.data.Sheet2(:,4)+MA2014.data.Sheet2(:,4)+DE2014.data.Sheet2(:,4)+NC2014.data.Sheet2(:,4)+CT2014.data.Sheet2(:,4)+RI2014.data.Sheet2(:,4)+VA2014.data.Sheet2(:,4)+NJ2014.data.Sheet2(:,4)+NY2014.data.Sheet2(:,4);
PredictedFish14L=MD2014.data.Sheet3(:,4)+MA2014.data.Sheet3(:,4)+DE2014.data.Sheet3(:,4)+NC2014.data.Sheet3(:,4)+CT2014.data.Sheet3(:,4)+RI2014.data.Sheet3(:,4)+VA2014.data.Sheet3(:,4)+NJ2014.data.Sheet3(:,4)+NY2014.data.Sheet3(:,4);
PredictedFish14C=MD2014.data.Sheet4(:,4)+MA2014.data.Sheet4(:,4)+DE2014.data.Sheet4(:,4)+NC2014.data.Sheet4(:,4)+CT2014.data.Sheet4(:,4)+RI2014.data.Sheet4(:,4)+VA2014.data.Sheet4(:,4)+NJ2014.data.Sheet4(:,4)+NY2014.data.Sheet4(:,4);

PredictedFish15A=MD2015.data.Sheet1(:,4)+MA2015.data.Sheet1(:,4)+DE2015.data.Sheet1(:,4)+NC2015.data.Sheet1(:,4)+CT2015.data.Sheet1(:,4)+RI2015.data.Sheet1(:,4)+VA2015.data.Sheet1(:,4)+NJ2015.data.Sheet1(:,4)+NY2015.data.Sheet1(:,4);
PredictedFish15M=MD2015.data.Sheet2(:,4)+MA2015.data.Sheet2(:,4)+DE2015.data.Sheet2(:,4)+NC2015.data.Sheet2(:,4)+CT2015.data.Sheet2(:,4)+RI2015.data.Sheet2(:,4)+VA2015.data.Sheet2(:,4)+NJ2015.data.Sheet2(:,4)+NY2015.data.Sheet2(:,4);
PredictedFish15L=MD2015.data.Sheet3(:,4)+MA2015.data.Sheet3(:,4)+DE2015.data.Sheet3(:,4)+NC2015.data.Sheet3(:,4)+CT2015.data.Sheet3(:,4)+RI2015.data.Sheet3(:,4)+VA2015.data.Sheet3(:,4)+NJ2015.data.Sheet3(:,4)+NY2015.data.Sheet3(:,4);
PredictedFish15C=MD2015.data.Sheet4(:,4)+MA2015.data.Sheet4(:,4)+DE2015.data.Sheet4(:,4)+NC2015.data.Sheet4(:,4)+CT2015.data.Sheet4(:,4)+RI2015.data.Sheet4(:,4)+VA2015.data.Sheet4(:,4)+NJ2015.data.Sheet4(:,4)+NY2015.data.Sheet4(:,4);

PredictedFish16A=MD2016.data.Sheet1(:,4)+MA2016.data.Sheet1(:,4)+DE2016.data.Sheet1(:,4)+NC2016.data.Sheet1(:,4)+CT2016.data.Sheet1(:,4)+RI2016.data.Sheet1(:,4)+VA2016.data.Sheet1(:,4)+NJ2016.data.Sheet1(:,4)+NY2016.data.Sheet1(:,4);
PredictedFish16M=MD2016.data.Sheet2(:,4)+MA2016.data.Sheet2(:,4)+DE2016.data.Sheet2(:,4)+NC2016.data.Sheet2(:,4)+CT2016.data.Sheet2(:,4)+RI2016.data.Sheet2(:,4)+VA2016.data.Sheet2(:,4)+NJ2016.data.Sheet2(:,4)+NY2016.data.Sheet2(:,4);
PredictedFish16L=MD2016.data.Sheet3(:,4)+MA2016.data.Sheet3(:,4)+DE2016.data.Sheet3(:,4)+NC2016.data.Sheet3(:,4)+CT2016.data.Sheet3(:,4)+RI2016.data.Sheet3(:,4)+VA2016.data.Sheet3(:,4)+NJ2016.data.Sheet3(:,4)+NY2016.data.Sheet3(:,4);
PredictedFish16C=MD2016.data.Sheet4(:,4)+MA2016.data.Sheet4(:,4)+DE2016.data.Sheet4(:,4)+NC2016.data.Sheet4(:,4)+CT2016.data.Sheet4(:,4)+RI2016.data.Sheet4(:,4)+VA2016.data.Sheet4(:,4)+NJ2016.data.Sheet4(:,4)+NY2016.data.Sheet4(:,4);

PredictedFish17A=MD2017.data.Sheet1(:,4)+MA2017.data.Sheet1(:,4)+DE2017.data.Sheet1(:,4)+NC2017.data.Sheet1(:,4)+CT2017.data.Sheet1(:,4)+RI2017.data.Sheet1(:,4)+VA2017.data.Sheet1(:,4)+NJ2017.data.Sheet1(:,4)+NY2017.data.Sheet1(:,4);
PredictedFish17M=MD2017.data.Sheet2(:,4)+MA2017.data.Sheet2(:,4)+DE2017.data.Sheet2(:,4)+NC2017.data.Sheet2(:,4)+CT2017.data.Sheet2(:,4)+RI2017.data.Sheet2(:,4)+VA2017.data.Sheet2(:,4)+NJ2017.data.Sheet2(:,4)+NY2017.data.Sheet2(:,4);
PredictedFish17L=MD2017.data.Sheet3(:,4)+MA2017.data.Sheet3(:,4)+DE2017.data.Sheet3(:,4)+NC2017.data.Sheet3(:,4)+CT2017.data.Sheet3(:,4)+RI2017.data.Sheet3(:,4)+VA2017.data.Sheet3(:,4)+NJ2017.data.Sheet3(:,4)+NY2017.data.Sheet3(:,4);
PredictedFish17C=MD2017.data.Sheet4(:,4)+MA2017.data.Sheet4(:,4)+DE2017.data.Sheet4(:,4)+NC2017.data.Sheet4(:,4)+CT2017.data.Sheet4(:,4)+RI2017.data.Sheet4(:,4)+VA2017.data.Sheet4(:,4)+NJ2017.data.Sheet4(:,4)+NY2017.data.Sheet4(:,4);

PredictedFish18A=MD2018.data.Sheet1(:,4)+MA2018.data.Sheet1(:,4)+DE2018.data.Sheet1(:,4)+NC2018.data.Sheet1(:,4)+CT2018.data.Sheet1(:,4)+RI2018.data.Sheet1(:,4)+VA2018.data.Sheet1(:,4)+NJ2018.data.Sheet1(:,4)+NY2018.data.Sheet1(:,4);
PredictedFish18M=MD2018.data.Sheet2(:,4)+MA2018.data.Sheet2(:,4)+DE2018.data.Sheet2(:,4)+NC2018.data.Sheet2(:,4)+CT2018.data.Sheet2(:,4)+RI2018.data.Sheet2(:,4)+VA2018.data.Sheet2(:,4)+NJ2018.data.Sheet2(:,4)+NY2018.data.Sheet2(:,4);
PredictedFish18L=MD2018.data.Sheet3(:,4)+MA2018.data.Sheet3(:,4)+DE2018.data.Sheet3(:,4)+NC2018.data.Sheet3(:,4)+CT2018.data.Sheet3(:,4)+RI2018.data.Sheet3(:,4)+VA2018.data.Sheet3(:,4)+NJ2018.data.Sheet3(:,4)+NY2018.data.Sheet3(:,4);
PredictedFish18C=MD2018.data.Sheet4(:,4)+MA2018.data.Sheet4(:,4)+DE2018.data.Sheet4(:,4)+NC2018.data.Sheet4(:,4)+CT2018.data.Sheet4(:,4)+RI2018.data.Sheet4(:,4)+VA2018.data.Sheet4(:,4)+NJ2018.data.Sheet4(:,4)+NY2018.data.Sheet4(:,4);

PredictedFish19A=MD2019.data.Sheet1(:,4)+MA2019.data.Sheet1(:,4)+DE2019.data.Sheet1(:,4)+NC2019.data.Sheet1(:,4)+CT2019.data.Sheet1(:,4)+RI2019.data.Sheet1(:,4)+VA2019.data.Sheet1(:,4)+NJ2019.data.Sheet1(:,4)+NY2019.data.Sheet1(:,4);
PredictedFish19M=MD2019.data.Sheet2(:,4)+MA2019.data.Sheet2(:,4)+DE2019.data.Sheet2(:,4)+NC2019.data.Sheet2(:,4)+CT2019.data.Sheet2(:,4)+RI2019.data.Sheet2(:,4)+VA2019.data.Sheet2(:,4)+NJ2019.data.Sheet2(:,4)+NY2019.data.Sheet2(:,4);
PredictedFish19L=MD2019.data.Sheet3(:,4)+MA2019.data.Sheet3(:,4)+DE2019.data.Sheet3(:,4)+NC2019.data.Sheet3(:,4)+CT2019.data.Sheet3(:,4)+RI2019.data.Sheet3(:,4)+VA2019.data.Sheet3(:,4)+NJ2019.data.Sheet3(:,4)+NY2019.data.Sheet3(:,4);
PredictedFish19C=MD2019.data.Sheet4(:,4)+MA2019.data.Sheet4(:,4)+DE2019.data.Sheet4(:,4)+NC2019.data.Sheet4(:,4)+CT2019.data.Sheet4(:,4)+RI2019.data.Sheet4(:,4)+VA2019.data.Sheet4(:,4)+NJ2019.data.Sheet4(:,4)+NY2019.data.Sheet4(:,4);

PredictedFish20A=MD2020.data.Sheet1(:,4)+MA2020.data.Sheet1(:,4)+DE2020.data.Sheet1(:,4)+NC2020.data.Sheet1(:,4)+CT2020.data.Sheet1(:,4)+RI2020.data.Sheet1(:,4)+VA2020.data.Sheet1(:,4)+NJ2020.data.Sheet1(:,4)+NY2020.data.Sheet1(:,4);
PredictedFish20M=MD2020.data.Sheet2(:,4)+MA2020.data.Sheet2(:,4)+DE2020.data.Sheet2(:,4)+NC2020.data.Sheet2(:,4)+CT2020.data.Sheet2(:,4)+RI2020.data.Sheet2(:,4)+VA2020.data.Sheet2(:,4)+NJ2020.data.Sheet2(:,4)+NY2020.data.Sheet2(:,4);
PredictedFish20L=MD2020.data.Sheet3(:,4)+MA2020.data.Sheet3(:,4)+DE2020.data.Sheet3(:,4)+NC2020.data.Sheet3(:,4)+CT2020.data.Sheet3(:,4)+RI2020.data.Sheet3(:,4)+VA2020.data.Sheet3(:,4)+NJ2020.data.Sheet3(:,4)+NY2020.data.Sheet3(:,4);
PredictedFish20C=MD2020.data.Sheet4(:,4)+MA2020.data.Sheet4(:,4)+DE2020.data.Sheet4(:,4)+NC2020.data.Sheet4(:,4)+CT2020.data.Sheet4(:,4)+RI2020.data.Sheet4(:,4)+VA2020.data.Sheet4(:,4)+NJ2020.data.Sheet4(:,4)+NY2020.data.Sheet4(:,4);


%Deviations on trips
predictedSEL=((PredictedTrips10L-ActualTrips10/1000).^2+(PredictedTrips11L-ActualTrips11/1000).^2+(PredictedTrips12L-ActualTrips12/1000).^2+(PredictedTrips13L-ActualTrips13/1000).^2+(PredictedTrips14L-ActualTrips14/1000).^2+(PredictedTrips15L-ActualTrips15/1000).^2+(PredictedTrips16L-ActualTrips16/1000).^2+(PredictedTrips17L-ActualTrips17/1000).^2+(PredictedTrips18L-ActualTrips18/1000).^2+(PredictedTrips19L-ActualTrips19/1000).^2+(PredictedTrips20L-ActualTrips20/1000).^2)/11;
predictedSEC=((PredictedTrips10C-ActualTrips10/1000).^2+(PredictedTrips11C-ActualTrips11/1000).^2+(PredictedTrips12C-ActualTrips12/1000).^2+(PredictedTrips13C-ActualTrips13/1000).^2+(PredictedTrips14C-ActualTrips14/1000).^2+(PredictedTrips15C-ActualTrips15/1000).^2+(PredictedTrips16C-ActualTrips16/1000).^2+(PredictedTrips17C-ActualTrips17/1000).^2+(PredictedTrips18C-ActualTrips18/1000).^2+(PredictedTrips19C-ActualTrips19/1000).^2+(PredictedTrips20C-ActualTrips20/1000).^2)/11;
predictedSEA=((PredictedTrips10A-ActualTrips10/1000).^2+(PredictedTrips11A-ActualTrips11/1000).^2+(PredictedTrips12A-ActualTrips12/1000).^2+(PredictedTrips13A-ActualTrips13/1000).^2+(PredictedTrips14A-ActualTrips14/1000).^2+(PredictedTrips15A-ActualTrips15/1000).^2+(PredictedTrips16A-ActualTrips16/1000).^2+(PredictedTrips17A-ActualTrips17/1000).^2+(PredictedTrips18A-ActualTrips18/1000).^2+(PredictedTrips19A-ActualTrips19/1000).^2+(PredictedTrips20A-ActualTrips20/1000).^2)/11;
predictedSEM=((PredictedTrips10M-ActualTrips10/1000).^2+(PredictedTrips11M-ActualTrips11/1000).^2+(PredictedTrips12M-ActualTrips12/1000).^2+(PredictedTrips13M-ActualTrips13/1000).^2+(PredictedTrips14M-ActualTrips14/1000).^2+(PredictedTrips15M-ActualTrips15/1000).^2+(PredictedTrips16M-ActualTrips16/1000).^2+(PredictedTrips17M-ActualTrips17/1000).^2+(PredictedTrips18M-ActualTrips18/1000).^2+(PredictedTrips19M-ActualTrips19/1000).^2+(PredictedTrips20M-ActualTrips20/1000).^2)/11;

MSEL=((mean(PredictedTrips10L)-ActualTrips10/1000)^2+(mean(PredictedTrips11L)-ActualTrips11/1000)^2+(mean(PredictedTrips12L)-ActualTrips12/1000)^2+(mean(PredictedTrips13L)-ActualTrips13/1000)^2+(mean(PredictedTrips14L)-ActualTrips14/1000)^2+(mean(PredictedTrips15L)-ActualTrips15/1000)^2+(mean(PredictedTrips16L)-ActualTrips16/1000)^2+(mean(PredictedTrips17L)-ActualTrips17/1000)^2+(mean(PredictedTrips18L)-ActualTrips18/1000)^2+(mean(PredictedTrips19L)-ActualTrips19/1000)^2+(mean(PredictedTrips20L)-ActualTrips20/1000)^2)/11;
MSEC=((mean(PredictedTrips10C)-ActualTrips10/1000)^2+(mean(PredictedTrips11C)-ActualTrips11/1000)^2+(mean(PredictedTrips12C)-ActualTrips12/1000)^2+(mean(PredictedTrips13C)-ActualTrips13/1000)^2+(mean(PredictedTrips14C)-ActualTrips14/1000)^2+(mean(PredictedTrips15C)-ActualTrips15/1000)^2+(mean(PredictedTrips16C)-ActualTrips16/1000)^2+(mean(PredictedTrips17C)-ActualTrips17/1000)^2+(mean(PredictedTrips18C)-ActualTrips18/1000)^2+(mean(PredictedTrips19C)-ActualTrips19/1000)^2+(mean(PredictedTrips20C)-ActualTrips20/1000)^2)/11;
MSEA=((mean(PredictedTrips10A)-ActualTrips10/1000)^2+(mean(PredictedTrips11A)-ActualTrips11/1000)^2+(mean(PredictedTrips12A)-ActualTrips12/1000)^2+(mean(PredictedTrips13A)-ActualTrips13/1000)^2+(mean(PredictedTrips14A)-ActualTrips14/1000)^2+(mean(PredictedTrips15A)-ActualTrips15/1000)^2+(mean(PredictedTrips16A)-ActualTrips16/1000)^2+(mean(PredictedTrips17A)-ActualTrips17/1000)^2+(mean(PredictedTrips18A)-ActualTrips18/1000)^2+(mean(PredictedTrips19A)-ActualTrips19/1000)^2+(mean(PredictedTrips20A)-ActualTrips20/1000)^2)/11;
MSEM=((mean(PredictedTrips10M)-ActualTrips10/1000)^2+(mean(PredictedTrips11M)-ActualTrips11/1000)^2+(mean(PredictedTrips12M)-ActualTrips12/1000)^2+(mean(PredictedTrips13M)-ActualTrips13/1000)^2+(mean(PredictedTrips14M)-ActualTrips14/1000)^2+(mean(PredictedTrips15M)-ActualTrips15/1000)^2+(mean(PredictedTrips16M)-ActualTrips16/1000)^2+(mean(PredictedTrips17M)-ActualTrips17/1000)^2+(mean(PredictedTrips18M)-ActualTrips18/1000)^2+(mean(PredictedTrips19M)-ActualTrips19/1000)^2+(mean(PredictedTrips20M)-ActualTrips20/1000)^2)/11;

%Calculating average deviation for predicted trips as a percentage of actual number of trips
DevSEL=100*horzcat(mean(abs((PredictedTrips10L-ActualTrips10/1000)))/(ActualTrips10/1000),mean(abs((PredictedTrips11L-ActualTrips11/1000)))/(ActualTrips11/1000), mean(abs((PredictedTrips12L-ActualTrips12/1000)))/(ActualTrips12/1000), mean(abs((PredictedTrips13L-ActualTrips13/1000)))/(ActualTrips13/1000), mean(abs((PredictedTrips14L-ActualTrips14/1000)))/(ActualTrips14/1000), mean(abs((PredictedTrips15L-ActualTrips15/1000)))/(ActualTrips15/1000), mean(abs((PredictedTrips16L-ActualTrips16/1000)))/(ActualTrips16/1000), mean(abs((PredictedTrips17L-ActualTrips17/1000)))/(ActualTrips17/1000), mean(abs((PredictedTrips18L-ActualTrips18/1000)))/(ActualTrips18/1000), mean(abs((PredictedTrips19L-ActualTrips19/1000)))/(ActualTrips19/1000), mean(abs((PredictedTrips20L-ActualTrips20/1000)))/(ActualTrips20/1000));
range_DevSEL=[min(DevSEL), max(DevSEL)];
meanDevSEL=mean(DevSEL);
DevSEC=100*horzcat(mean(abs((PredictedTrips10C-ActualTrips10/1000)))/(ActualTrips10/1000),mean(abs((PredictedTrips11C-ActualTrips11/1000)))/(ActualTrips11/1000), mean(abs((PredictedTrips12C-ActualTrips12/1000)))/(ActualTrips12/1000), mean(abs((PredictedTrips13C-ActualTrips13/1000)))/(ActualTrips13/1000), mean(abs((PredictedTrips14C-ActualTrips14/1000)))/(ActualTrips14/1000), mean(abs((PredictedTrips15C-ActualTrips15/1000)))/(ActualTrips15/1000), mean(abs((PredictedTrips16C-ActualTrips16/1000)))/(ActualTrips16/1000), mean(abs((PredictedTrips17C-ActualTrips17/1000)))/(ActualTrips17/1000), mean(abs((PredictedTrips18C-ActualTrips18/1000)))/(ActualTrips18/1000), mean(abs((PredictedTrips19C-ActualTrips19/1000)))/(ActualTrips19/1000), mean(abs((PredictedTrips20C-ActualTrips20/1000)))/(ActualTrips20/1000));
range_DevSEC=[min(DevSEC), max(DevSEC)];
meanDevSEC=mean(DevSEC);
DevSEA=100*horzcat(mean(abs((PredictedTrips10A-ActualTrips10/1000)))/(ActualTrips10/1000),mean(abs((PredictedTrips11A-ActualTrips11/1000)))/(ActualTrips11/1000), mean(abs((PredictedTrips12A-ActualTrips12/1000)))/(ActualTrips12/1000), mean(abs((PredictedTrips13A-ActualTrips13/1000)))/(ActualTrips13/1000), mean(abs((PredictedTrips14A-ActualTrips14/1000)))/(ActualTrips14/1000), mean(abs((PredictedTrips15A-ActualTrips15/1000)))/(ActualTrips15/1000), mean(abs((PredictedTrips16A-ActualTrips16/1000)))/(ActualTrips16/1000), mean(abs((PredictedTrips17A-ActualTrips17/1000)))/(ActualTrips17/1000), mean(abs((PredictedTrips18A-ActualTrips18/1000)))/(ActualTrips18/1000), mean(abs((PredictedTrips19A-ActualTrips19/1000)))/(ActualTrips19/1000), mean(abs((PredictedTrips20A-ActualTrips20/1000)))/(ActualTrips20/1000));
range_DevSEA=[min(DevSEA), max(DevSEA)];
meanDevSEA=mean(DevSEA);
DevSEM=100*horzcat(mean(abs((PredictedTrips10M-ActualTrips10/1000)))/(ActualTrips10/1000),mean(abs((PredictedTrips11M-ActualTrips11/1000)))/(ActualTrips11/1000), mean(abs((PredictedTrips12M-ActualTrips12/1000)))/(ActualTrips12/1000), mean(abs((PredictedTrips13M-ActualTrips13/1000)))/(ActualTrips13/1000), mean(abs((PredictedTrips14M-ActualTrips14/1000)))/(ActualTrips14/1000), mean(abs((PredictedTrips15M-ActualTrips15/1000)))/(ActualTrips15/1000), mean(abs((PredictedTrips16M-ActualTrips16/1000)))/(ActualTrips16/1000), mean(abs((PredictedTrips17M-ActualTrips17/1000)))/(ActualTrips17/1000), mean(abs((PredictedTrips18M-ActualTrips18/1000)))/(ActualTrips18/1000), mean(abs((PredictedTrips19M-ActualTrips19/1000)))/(ActualTrips19/1000), mean(abs((PredictedTrips20M-ActualTrips20/1000)))/(ActualTrips20/1000));
range_DevSEM=[min(DevSEM), max(DevSEM)];
meanDevSEM=mean(DevSEM);

%Deviations on landings
predictedSELFish=((PredictedFish10L-ActualFish10/1000).^2+(PredictedFish11L-ActualFish11/1000).^2+(PredictedFish12L-ActualFish12/1000).^2+(PredictedFish13L-ActualFish13/1000).^2+(PredictedFish14L-ActualFish14/1000).^2+(PredictedFish15L-ActualFish15/1000).^2+(PredictedFish16L-ActualFish16/1000).^2+(PredictedFish17L-ActualFish17/1000).^2+(PredictedFish18L-ActualFish18/1000).^2+(PredictedFish19L-ActualFish19/1000).^2+(PredictedFish20L-ActualFish20/1000).^2)/11;
predictedSECFish=((PredictedFish10C-ActualFish10/1000).^2+(PredictedFish11C-ActualFish11/1000).^2+(PredictedFish12C-ActualFish12/1000).^2+(PredictedFish13C-ActualFish13/1000).^2+(PredictedFish14C-ActualFish14/1000).^2+(PredictedFish15C-ActualFish15/1000).^2+(PredictedFish16C-ActualFish16/1000).^2+(PredictedFish17C-ActualFish17/1000).^2+(PredictedFish18C-ActualFish18/1000).^2+(PredictedFish19C-ActualFish19/1000).^2+(PredictedFish20C-ActualFish20/1000).^2)/11;
predictedSEAFish=((PredictedFish10A-ActualFish10/1000).^2+(PredictedFish11A-ActualFish11/1000).^2+(PredictedFish12A-ActualFish12/1000).^2+(PredictedFish13A-ActualFish13/1000).^2+(PredictedFish14A-ActualFish14/1000).^2+(PredictedFish15A-ActualFish15/1000).^2+(PredictedFish16A-ActualFish16/1000).^2+(PredictedFish17A-ActualFish17/1000).^2+(PredictedFish18A-ActualFish18/1000).^2+(PredictedFish19A-ActualFish19/1000).^2+(PredictedFish20A-ActualFish20/1000).^2)/11;
predictedSEMFish=((PredictedFish10M-ActualFish10/1000).^2+(PredictedFish11M-ActualFish11/1000).^2+(PredictedFish12M-ActualFish12/1000).^2+(PredictedFish13M-ActualFish13/1000).^2+(PredictedFish14M-ActualFish14/1000).^2+(PredictedFish15M-ActualFish15/1000).^2+(PredictedFish16M-ActualFish16/1000).^2+(PredictedFish17M-ActualFish17/1000).^2+(PredictedFish18M-ActualFish18/1000).^2+(PredictedFish19M-ActualFish19/1000).^2+(PredictedFish20M-ActualFish20/1000).^2)/11;

MSELFish=((mean(PredictedFish10L)-ActualFish10/1000)^2+(mean(PredictedFish11L)-ActualFish11/1000)^2+(mean(PredictedFish12L)-ActualFish12/1000)^2+(mean(PredictedFish13L)-ActualFish13/1000)^2+(mean(PredictedFish14L)-ActualFish14/1000)^2+(mean(PredictedFish15L)-ActualFish15/1000)^2+(mean(PredictedFish16L)-ActualFish16/1000)^2+(mean(PredictedFish17L)-ActualFish17/1000)^2+(mean(PredictedFish18L)-ActualFish18/1000)^2+(mean(PredictedFish19L)-ActualFish19/1000)^2+(mean(PredictedFish20L)-ActualFish20/1000)^2)/11;
MSECFish=((mean(PredictedFish10C)-ActualFish10/1000)^2+(mean(PredictedFish11C)-ActualFish11/1000)^2+(mean(PredictedFish12C)-ActualFish12/1000)^2+(mean(PredictedFish13C)-ActualFish13/1000)^2+(mean(PredictedFish14C)-ActualFish14/1000)^2+(mean(PredictedFish15C)-ActualFish15/1000)^2+(mean(PredictedFish16C)-ActualFish16/1000)^2+(mean(PredictedFish17C)-ActualFish17/1000)^2+(mean(PredictedFish18C)-ActualFish18/1000)^2+(mean(PredictedFish19C)-ActualFish19/1000)^2+(mean(PredictedFish20C)-ActualFish20/1000)^2)/11;
MSEAFish=((mean(PredictedFish10A)-ActualFish10/1000)^2+(mean(PredictedFish11A)-ActualFish11/1000)^2+(mean(PredictedFish12A)-ActualFish12/1000)^2+(mean(PredictedFish13A)-ActualFish13/1000)^2+(mean(PredictedFish14A)-ActualFish14/1000)^2+(mean(PredictedFish15A)-ActualFish15/1000)^2+(mean(PredictedFish16A)-ActualFish16/1000)^2+(mean(PredictedFish17A)-ActualFish17/1000)^2+(mean(PredictedFish18A)-ActualFish18/1000)^2+(mean(PredictedFish19A)-ActualFish19/1000)^2+(mean(PredictedFish20A)-ActualFish20/1000)^2)/11;
MSEMFish=((mean(PredictedFish10M)-ActualFish10/1000)^2+(mean(PredictedFish11M)-ActualFish11/1000)^2+(mean(PredictedFish12M)-ActualFish12/1000)^2+(mean(PredictedFish13M)-ActualFish13/1000)^2+(mean(PredictedFish14M)-ActualFish14/1000)^2+(mean(PredictedFish15M)-ActualFish15/1000)^2+(mean(PredictedFish16M)-ActualFish16/1000)^2+(mean(PredictedFish17M)-ActualFish17/1000)^2+(mean(PredictedFish18M)-ActualFish18/1000)^2+(mean(PredictedFish19M)-ActualFish19/1000)^2+(mean(PredictedFish20M)-ActualFish20/1000)^2)/11;

%Calculating average deviation for predicted landings a percentage of actual landings
DevSEFishL=100*horzcat(mean(abs((PredictedFish10L-ActualFish10/1000)))/(ActualFish10/1000),mean(abs((PredictedFish11L-ActualFish11/1000)))/(ActualFish11/1000), mean(abs((PredictedFish12L-ActualFish12/1000)))/(ActualFish12/1000), mean(abs((PredictedFish13L-ActualFish13/1000)))/(ActualFish13/1000), mean(abs((PredictedFish14L-ActualFish14/1000)))/(ActualFish14/1000), mean(abs((PredictedFish15L-ActualFish15/1000)))/(ActualFish15/1000), mean(abs((PredictedFish16L-ActualFish16/1000)))/(ActualFish16/1000), mean(abs((PredictedFish17L-ActualFish17/1000)))/(ActualFish17/1000), mean(abs((PredictedFish18L-ActualFish18/1000)))/(ActualFish18/1000), mean(abs((PredictedFish19L-ActualFish19/1000)))/(ActualFish19/1000), mean(abs((PredictedFish20L-ActualFish20/1000)))/(ActualFish20/1000));
range_DevSEFishL=[min(DevSEFishL), max(DevSEFishL)];
meanDevSEFishL=mean(DevSEFishL);
DevSEFishC=100*horzcat(mean(abs((PredictedFish10C-ActualFish10/1000)))/(ActualFish10/1000),mean(abs((PredictedFish11C-ActualFish11/1000)))/(ActualFish11/1000), mean(abs((PredictedFish12C-ActualFish12/1000)))/(ActualFish12/1000), mean(abs((PredictedFish13C-ActualFish13/1000)))/(ActualFish13/1000), mean(abs((PredictedFish14C-ActualFish14/1000)))/(ActualFish14/1000), mean(abs((PredictedFish15C-ActualFish15/1000)))/(ActualFish15/1000), mean(abs((PredictedFish16C-ActualFish16/1000)))/(ActualFish16/1000), mean(abs((PredictedFish17C-ActualFish17/1000)))/(ActualFish17/1000), mean(abs((PredictedFish18C-ActualFish18/1000)))/(ActualFish18/1000), mean(abs((PredictedFish19C-ActualFish19/1000)))/(ActualFish19/1000), mean(abs((PredictedFish20C-ActualFish20/1000)))/(ActualFish20/1000));
range_DevSEFishC=[min(DevSEFishC), max(DevSEFishC)];
meanDevSEFishC=mean(DevSEFishC);
DevSEFishA=100*horzcat(mean(abs((PredictedFish10A-ActualFish10/1000)))/(ActualFish10/1000),mean(abs((PredictedFish11A-ActualFish11/1000)))/(ActualFish11/1000), mean(abs((PredictedFish12A-ActualFish12/1000)))/(ActualFish12/1000), mean(abs((PredictedFish13A-ActualFish13/1000)))/(ActualFish13/1000), mean(abs((PredictedFish14A-ActualFish14/1000)))/(ActualFish14/1000), mean(abs((PredictedFish15A-ActualFish15/1000)))/(ActualFish15/1000), mean(abs((PredictedFish16A-ActualFish16/1000)))/(ActualFish16/1000), mean(abs((PredictedFish17A-ActualFish17/1000)))/(ActualFish17/1000), mean(abs((PredictedFish18A-ActualFish18/1000)))/(ActualFish18/1000), mean(abs((PredictedFish19A-ActualFish19/1000)))/(ActualFish19/1000), mean(abs((PredictedFish20A-ActualFish20/1000)))/(ActualFish20/1000));
range_DevSEFishA=[min(DevSEFishA), max(DevSEFishA)];
meanDevSEFishA=mean(DevSEFishA);
DevSEFishM=100*horzcat(mean(abs((PredictedFish10M-ActualFish10/1000)))/(ActualFish10/1000),mean(abs((PredictedFish11M-ActualFish11/1000)))/(ActualFish11/1000), mean(abs((PredictedFish12M-ActualFish12/1000)))/(ActualFish12/1000), mean(abs((PredictedFish13M-ActualFish13/1000)))/(ActualFish13/1000), mean(abs((PredictedFish14M-ActualFish14/1000)))/(ActualFish14/1000), mean(abs((PredictedFish15M-ActualFish15/1000)))/(ActualFish15/1000), mean(abs((PredictedFish16M-ActualFish16/1000)))/(ActualFish16/1000), mean(abs((PredictedFish17M-ActualFish17/1000)))/(ActualFish17/1000), mean(abs((PredictedFish18M-ActualFish18/1000)))/(ActualFish18/1000), mean(abs((PredictedFish19M-ActualFish19/1000)))/(ActualFish19/1000), mean(abs((PredictedFish20M-ActualFish20/1000)))/(ActualFish20/1000));
range_DevSEFishM=[min(DevSEFishM), max(DevSEFishM)];
meanDevSEFishM=mean(DevSEFishM);

%Figure 1 and Table 3
figure
subplot(2,1,1)
h(1,1)=cdfplot(predictedSEA);
hold on 
h(1,2)=cdfplot(predictedSEM);
hold on 
h(1,3)=cdfplot(predictedSEL);
hold on 
h(1,4)=cdfplot(predictedSEC);
title('(a)','FontSize',14)
legend('$\alpha$-MEU','mLogit','clogit','CARA','Location','best','Interpreter','latex')
xlabel('MSPE (Trips)', 'FontSize',18)
set(gca,'FontSize',18)
set(h(:,1), 'LineStyle', '-', 'Color', 'k','LineWidth',2);
set(h(:,2), 'LineStyle', '-.', 'Color', [0.4 0.4 0.4], 'LineWidth',2);
set(h(:,3), 'LineStyle', '--', 'Color', 'k','LineWidth',2);
set(h(:,4), 'LineStyle', ':', 'Color', [0.6 0.6 0.6], 'LineWidth',2);
hold off

subplot(2,1,2)
h(1,1)=cdfplot(predictedSEAFish);
hold on 
h(1,2)=cdfplot(predictedSEMFish);
hold on 
h(1,3)=cdfplot(predictedSELFish);
hold on 
h(1,4)=cdfplot(predictedSECFish);
title('(b)','FontSize',14)
legend('$\alpha$-MEU','mLogit','clogit','CARA','Location','best','Interpreter','latex')
xlabel('MSPE (Landings)', 'FontSize',18)
set(gca,'FontSize',18)
set(h(:,1), 'LineStyle', '-', 'Color', 'k','LineWidth',2);
set(h(:,2), 'LineStyle', '-.', 'Color', [0.4 0.4 0.4], 'LineWidth',2);
set(h(:,3), 'LineStyle', '--', 'Color', 'k','LineWidth',2);
set(h(:,4), 'LineStyle', ':', 'Color', [0.6 0.6 0.6], 'LineWidth',2);
hold off
f=gcf;
savefig(gcf,'Figure1.fig');

%FSD tests
%Trips
pval_trips1=TestFSD(predictedSEM,predictedSEA,'Pooled');
pval_trips2=TestFSD(predictedSEL,predictedSEA,'Pooled');
pval_trips3=TestFSD(predictedSEC,predictedSEA,'Pooled');
%Landings
pval_landings1=TestFSD(predictedSEMFish,predictedSEAFish,'Pooled');
pval_landings2=TestFSD(predictedSELFish,predictedSEAFish,'Pooled');
pval_landings3=TestFSD(predictedSECFish,predictedSEAFish,'Pooled');

Table3=vertcat(horzcat("Null Hypothesis (H0)","p-value (1,000 bootstraps)", "FSD Observed"),horzcat("MSPE (Trips)","",""),vertcat(horzcat("mlogit FSD alpha-maxmin",pval_trips1,pval_trips1>0.05), horzcat("clogit FSD alpha-maxmin",pval_trips2,pval_trips2>0.05),horzcat("CARA FSD alpha-maxmin",pval_trips3, pval_trips3>0.05)),horzcat("MSPE (Landings)","",""),vertcat(horzcat("mlogit FSD alpha-maxmin",pval_landings1,pval_trips1>0.05),horzcat("clogit FSD alpha-maxmin",pval_landings2,pval_landings2>0.05),horzcat("CARA FSD alpha-maxmin",pval_landings3,pval_landings3>0.05)));

%Two-sample t-tests for equal group means in Table A-8
%Trips
[h20a,p20a]=ttest2(PredictedTrips20M,PredictedTrips20A);
[h20b,p20b]=ttest2(PredictedTrips20C,PredictedTrips20A);
[h20c,p20c]=ttest2(PredictedTrips20L,PredictedTrips20A);

[h19a,p19a]=ttest2(PredictedTrips19M,PredictedTrips19A);
[h19b,p19b]=ttest2(PredictedTrips19C,PredictedTrips19A);
[h19c,p19c]=ttest2(PredictedTrips19L,PredictedTrips19A);

[h18a,p18a]=ttest2(PredictedTrips18M,PredictedTrips18A);
[h18b,p18b]=ttest2(PredictedTrips18C,PredictedTrips18A);
[h18c,p18c]=ttest2(PredictedTrips18L,PredictedTrips18A);

[h17a,p17a]=ttest2(PredictedTrips17M,PredictedTrips17A);
[h17b,p17b]=ttest2(PredictedTrips17C,PredictedTrips17A);
[h17c,p17c]=ttest2(PredictedTrips17L,PredictedTrips17A);

[h16a,p16a]=ttest2(PredictedTrips16M,PredictedTrips16A);
[h16b,p16b]=ttest2(PredictedTrips16C,PredictedTrips16A);
[h16c,p16c]=ttest2(PredictedTrips16L,PredictedTrips16A);

[h15a,p15a]=ttest2(PredictedTrips15M,PredictedTrips15A);
[h15b,p15b]=ttest2(PredictedTrips15C,PredictedTrips15A);
[h15c,p15c]=ttest2(PredictedTrips15L,PredictedTrips15A);

[h14a,p14a]=ttest2(PredictedTrips14M,PredictedTrips14A);
[h14b,p14b]=ttest2(PredictedTrips14C,PredictedTrips14A);
[h14c,p14c]=ttest2(PredictedTrips14L,PredictedTrips14A);

[h13a,p13a]=ttest2(PredictedTrips13M,PredictedTrips13A);
[h13b,p13b]=ttest2(PredictedTrips13C,PredictedTrips13A);
[h13c,p13c]=ttest2(PredictedTrips13L,PredictedTrips13A);

[h12a,p12a]=ttest2(PredictedTrips12M,PredictedTrips12A);
[h12b,p12b]=ttest2(PredictedTrips12C,PredictedTrips12A);
[h12c,p12c]=ttest2(PredictedTrips12L,PredictedTrips12A);

[h11a,p11a]=ttest2(PredictedTrips11M,PredictedTrips11A);
[h11b,p11b]=ttest2(PredictedTrips11C,PredictedTrips11A);
[h11c,p11c]=ttest2(PredictedTrips11L,PredictedTrips11A);

[h10a,p10a]=ttest2(PredictedTrips10M,PredictedTrips10A);
[h10b,p10b]=ttest2(PredictedTrips10C,PredictedTrips10A);
[h10c,p10c]=ttest2(PredictedTrips10L,PredictedTrips10A);

%Landings
[h20Fa,p20Fa]=ttest2(PredictedFish20M,PredictedFish20A);
[h20Fb,p20Fb]=ttest2(PredictedFish20C,PredictedFish20A);
[h20Fc,p20Fc]=ttest2(PredictedFish20L,PredictedFish20A);

[h19Fa,p19Fa]=ttest2(PredictedFish19M,PredictedFish19A);
[h19Fb,p19Fb]=ttest2(PredictedFish19C,PredictedFish19A);
[h19Fc,p19Fc]=ttest2(PredictedFish19L,PredictedFish19A);

[h18Fa,p18Fa]=ttest2(PredictedFish18M,PredictedFish18A);
[h18Fb,p18Fb]=ttest2(PredictedFish18C,PredictedFish18A);
[h18Fc,p18Fc]=ttest2(PredictedFish18L,PredictedFish18A);

[h17Fa,p17Fa]=ttest2(PredictedFish17M,PredictedFish17A);
[h17Fb,p17Fb]=ttest2(PredictedFish17C,PredictedFish17A);
[h17Fc,p17Fc]=ttest2(PredictedFish17L,PredictedFish17A);

[h16Fa,p16Fa]=ttest2(PredictedFish16M,PredictedFish16A);
[h16Fb,p16Fb]=ttest2(PredictedFish16C,PredictedFish16A);
[h16Fc,p16Fc]=ttest2(PredictedFish16L,PredictedFish16A);

[h15Fa,p15Fa]=ttest2(PredictedFish15M,PredictedFish15A);
[h15Fb,p15Fb]=ttest2(PredictedFish15C,PredictedFish15A);
[h15Fc,p15Fc]=ttest2(PredictedFish15L,PredictedFish15A);

[h14Fa,p14Fa]=ttest2(PredictedFish14M,PredictedFish14A);
[h14Fb,p14Fb]=ttest2(PredictedFish14C,PredictedFish14A);
[h14Fc,p14Fc]=ttest2(PredictedFish14L,PredictedFish14A);

[h13Fa,p13Fa]=ttest2(PredictedFish13M,PredictedFish13A);
[h13Fb,p13Fb]=ttest2(PredictedFish13C,PredictedFish13A);
[h13Fc,p13Fc]=ttest2(PredictedFish13L,PredictedFish13A);

[h12Fa,p12Fa]=ttest2(PredictedFish12M,PredictedFish12A);
[h12Fb,p12Fb]=ttest2(PredictedFish12C,PredictedFish12A);
[h12Fc,p12Fc]=ttest2(PredictedFish12L,PredictedFish12A);

[h11Fa,p11Fa]=ttest2(PredictedFish11M,PredictedFish11A);
[h11Fb,p11Fb]=ttest2(PredictedFish11C,PredictedFish11A);
[h11Fc,p11Fc]=ttest2(PredictedFish11L,PredictedFish11A);

[h10Fa,p10Fa]=ttest2(PredictedFish10M,PredictedFish10A);
[h10Fb,p10Fb]=ttest2(PredictedFish10C,PredictedFish10A);
[h10Fc,p10Fc]=ttest2(PredictedFish10L,PredictedFish10A);

body=horzcat(vertcat("2010","2011","2012","2013","2014","2015","2016","2017","2018","2019","2020"),vertcat(horzcat(p10a,p10b,p10c),horzcat(p11a,p11b,p11c),horzcat(p12a,p12b,p12c),horzcat(p13a,p13b,p13c),horzcat(p14a,p14b,p14c),horzcat(p15a,p15b,p15c),horzcat(p16a,p16b,p16c),horzcat(p17a,p17b,p17c),horzcat(p18a,p18b,p18c),horzcat(p19a,p19b,p19c),horzcat(p20a,p20b,p20c)));
heading=horzcat("Linear(mlogit)","CARA","Linear(clogit)");
bodyF=horzcat(vertcat("2010","2011","2012","2013","2014","2015","2016","2017","2018","2019","2020"),vertcat(horzcat(p10Fa,p10Fb,p10Fc),horzcat(p11Fa,p11Fb,p11Fc),horzcat(p12Fa,p12Fb,p12Fc),horzcat(p13Fa,p13Fb,p13Fc),horzcat(p14Fa,p14Fb,p14Fc),horzcat(p15Fa,p15Fb,p15Fc),horzcat(p16Fa,p16Fb,p16Fc),horzcat(p17Fa,p17Fb,p17Fc),horzcat(p18Fa,p18Fb,p18Fc),horzcat(p19Fa,p19Fb,p19Fc),horzcat(p20Fa,p20Fb,p20Fc)));
TableA8=vertcat(vertcat(horzcat("Trips",heading),body),vertcat(horzcat("Landings",heading),bodyF));
clearvars body bodyF

%Preparing data for the Stata file violin_plots.do that generates the
%violin plots (Figures A-4 through A-7)
models=vertcat(repmat("clogit",length(PredictedTrips10L),1),repmat("CARA",length(PredictedTrips10C),1),repmat("ambiguity",length(PredictedTrips10A),1),repmat("mlogit",length(PredictedTrips10M),1));
predictionTrips10=vertcat(PredictedTrips10L,PredictedTrips10C,PredictedTrips10A,PredictedTrips10M);
predictionTrips11=vertcat(PredictedTrips11L,PredictedTrips11C,PredictedTrips11A,PredictedTrips11M);
predictionTrips12=vertcat(PredictedTrips12L,PredictedTrips12C,PredictedTrips12A,PredictedTrips12M);
predictionTrips13=vertcat(PredictedTrips13L,PredictedTrips13C,PredictedTrips13A,PredictedTrips13M);
predictionTrips14=vertcat(PredictedTrips14L,PredictedTrips14C,PredictedTrips14A,PredictedTrips14M);
predictionTrips15=vertcat(PredictedTrips15L,PredictedTrips15C,PredictedTrips15A,PredictedTrips15M);
predictionTrips16=vertcat(PredictedTrips16L,PredictedTrips16C,PredictedTrips16A,PredictedTrips16M);
predictionTrips17=vertcat(PredictedTrips17L,PredictedTrips17C,PredictedTrips17A,PredictedTrips17M);
predictionTrips18=vertcat(PredictedTrips18L,PredictedTrips18C,PredictedTrips18A,PredictedTrips18M);
predictionTrips19=vertcat(PredictedTrips19L,PredictedTrips19C,PredictedTrips19A,PredictedTrips19M);
predictionTrips20=vertcat(PredictedTrips20L,PredictedTrips20C,PredictedTrips20A,PredictedTrips20M);
actualT10=ActualTrips10/1000*ones(4*length(PredictedTrips10L),1);
actualT11=ActualTrips11/1000*ones(4*length(PredictedTrips11L),1);
actualT12=ActualTrips12/1000*ones(4*length(PredictedTrips12L),1);
actualT13=ActualTrips13/1000*ones(4*length(PredictedTrips13L),1);
actualT14=ActualTrips14/1000*ones(4*length(PredictedTrips14L),1);
actualT15=ActualTrips15/1000*ones(4*length(PredictedTrips15L),1);
actualT16=ActualTrips16/1000*ones(4*length(PredictedTrips16L),1);
actualT17=ActualTrips17/1000*ones(4*length(PredictedTrips17L),1);
actualT18=ActualTrips18/1000*ones(4*length(PredictedTrips18L),1);
actualT19=ActualTrips19/1000*ones(4*length(PredictedTrips19L),1);
actualT20=ActualTrips20/1000*ones(4*length(PredictedTrips20L),1);
predictionFish10=vertcat(PredictedFish10L,PredictedFish10C,PredictedFish10A,PredictedFish10M);
predictionFish11=vertcat(PredictedFish11L,PredictedFish11C,PredictedFish11A,PredictedFish11M);
predictionFish12=vertcat(PredictedFish12L,PredictedFish12C,PredictedFish12A,PredictedFish12M);
predictionFish13=vertcat(PredictedFish13L,PredictedFish13C,PredictedFish13A,PredictedFish13M);
predictionFish14=vertcat(PredictedFish14L,PredictedFish14C,PredictedFish14A,PredictedFish14M);
predictionFish15=vertcat(PredictedFish15L,PredictedFish15C,PredictedFish15A,PredictedFish15M);
predictionFish16=vertcat(PredictedFish16L,PredictedFish16C,PredictedFish16A,PredictedFish16M);
predictionFish17=vertcat(PredictedFish17L,PredictedFish17C,PredictedFish17A,PredictedFish17M);
predictionFish18=vertcat(PredictedFish18L,PredictedFish18C,PredictedFish18A,PredictedFish18M);
predictionFish19=vertcat(PredictedFish19L,PredictedFish19C,PredictedFish19A,PredictedFish19M);
predictionFish20=vertcat(PredictedFish20L,PredictedFish20C,PredictedFish20A,PredictedFish20M);
actualF10=ActualFish10/1000*ones(4*length(PredictedFish10L),1);
actualF11=ActualFish11/1000*ones(4*length(PredictedFish11L),1);
actualF12=ActualFish12/1000*ones(4*length(PredictedFish12L),1);
actualF13=ActualFish13/1000*ones(4*length(PredictedFish13L),1);
actualF14=ActualFish14/1000*ones(4*length(PredictedFish14L),1);
actualF15=ActualFish15/1000*ones(4*length(PredictedFish15L),1);
actualF16=ActualFish16/1000*ones(4*length(PredictedFish16L),1);
actualF17=ActualFish17/1000*ones(4*length(PredictedFish17L),1);
actualF18=ActualFish18/1000*ones(4*length(PredictedFish18L),1);
actualF19=ActualFish19/1000*ones(4*length(PredictedFish19L),1);
actualF20=ActualFish20/1000*ones(4*length(PredictedFish20L),1);

varNames1={'models','PredT10','ActualT10','PredT11','ActualT11','PredT12','ActualT12','PredT13','ActualT13','PredT14','ActualT14','PredT15','ActualT15','PredT16','ActualT16','PredT17','ActualT17','PredT18','ActualT18','PredT19', 'ActualT19','PredT20','ActualT20','PredF10','ActualF10','PredF11','ActualF11','PredF12','ActualF12','PredF13','ActualF13','PredF14','ActualF14','PredF15','ActualF15','PredF16','ActualF16','PredF17','ActualF17','PredF18','ActualF18','PredF19','ActualF19','PredF20','ActualF20',};
writetable(table(models,predictionTrips10,actualT10,predictionTrips11,actualT11,predictionTrips12,actualT12,predictionTrips13,actualT13,predictionTrips14,actualT14,predictionTrips15,actualT15,predictionTrips16,actualT16,predictionTrips17,actualT17,predictionTrips18,actualT18,predictionTrips19,actualT19,predictionTrips20,actualT20,predictionFish10,actualF10,predictionFish11,actualF11,predictionFish12,actualF12,predictionFish13,actualF13,predictionFish14,actualF14,predictionFish15,actualF15,predictionFish16,actualF16,predictionFish17,actualF17,predictionFish18,actualF18,predictionFish19,actualF19,predictionFish20,actualF20,'VariableNames',varNames1),'G:\My Drive\RESEARCH\WORKING PAPERS\McConnell_Holzer on Responses to Risky Catch\Matlab code\simulations\Figures_forward\dataGraphs.xls','Sheet',1)

%Table 4 (MSPE for each model as % of MSPE for the preferred model)
% Trips
MSPEA_CT=((CT2010.data.Sheet1(:,1)-TripsCT10/1000).^2+(CT2011.data.Sheet1(:,1)-TripsCT11/1000).^2+(CT2012.data.Sheet1(:,1)-TripsCT12/1000).^2+(CT2013.data.Sheet1(:,1)-TripsCT13/1000).^2+(CT2014.data.Sheet1(:,1)-TripsCT14/1000).^2+(CT2015.data.Sheet1(:,1)-TripsCT15/1000).^2+(CT2016.data.Sheet1(:,1)-TripsCT16/1000).^2+(CT2017.data.Sheet1(:,1)-TripsCT17/1000).^2+(CT2018.data.Sheet1(:,1)-TripsCT18/1000).^2+(CT2019.data.Sheet1(:,1)-TripsCT19/1000).^2+(CT2020.data.Sheet1(:,1)-TripsCT20/1000).^2)/11;
MSPEM_CT=((CT2010.data.Sheet2(:,1)-TripsCT10/1000).^2+(CT2011.data.Sheet2(:,1)-TripsCT11/1000).^2+(CT2012.data.Sheet2(:,1)-TripsCT12/1000).^2+(CT2013.data.Sheet2(:,1)-TripsCT13/1000).^2+(CT2014.data.Sheet2(:,1)-TripsCT14/1000).^2+(CT2015.data.Sheet2(:,1)-TripsCT15/1000).^2+(CT2016.data.Sheet2(:,1)-TripsCT16/1000).^2+(CT2017.data.Sheet2(:,1)-TripsCT17/1000).^2+(CT2018.data.Sheet2(:,1)-TripsCT18/1000).^2+(CT2019.data.Sheet2(:,1)-TripsCT19/1000).^2+(CT2020.data.Sheet2(:,1)-TripsCT20/1000).^2)/11;
MSPEL_CT=((CT2010.data.Sheet3(:,1)-TripsCT10/1000).^2+(CT2011.data.Sheet3(:,1)-TripsCT11/1000).^2+(CT2012.data.Sheet3(:,1)-TripsCT12/1000).^2+(CT2013.data.Sheet3(:,1)-TripsCT13/1000).^2+(CT2014.data.Sheet3(:,1)-TripsCT14/1000).^2+(CT2015.data.Sheet3(:,1)-TripsCT15/1000).^2+(CT2016.data.Sheet3(:,1)-TripsCT16/1000).^2+(CT2017.data.Sheet3(:,1)-TripsCT17/1000).^2+(CT2018.data.Sheet3(:,1)-TripsCT18/1000).^2+(CT2019.data.Sheet3(:,1)-TripsCT19/1000).^2+(CT2020.data.Sheet3(:,1)-TripsCT20/1000).^2)/11;
MSPEC_CT=((CT2010.data.Sheet4(:,1)-TripsCT10/1000).^2+(CT2011.data.Sheet4(:,1)-TripsCT11/1000).^2+(CT2012.data.Sheet4(:,1)-TripsCT12/1000).^2+(CT2013.data.Sheet4(:,1)-TripsCT13/1000).^2+(CT2014.data.Sheet4(:,1)-TripsCT14/1000).^2+(CT2015.data.Sheet4(:,1)-TripsCT15/1000).^2+(CT2016.data.Sheet4(:,1)-TripsCT16/1000).^2+(CT2017.data.Sheet4(:,1)-TripsCT17/1000).^2+(CT2018.data.Sheet4(:,1)-TripsCT18/1000).^2+(CT2019.data.Sheet4(:,1)-TripsCT19/1000).^2+(CT2020.data.Sheet4(:,1)-TripsCT20/1000).^2)/11;

MSPEA_DE=((DE2010.data.Sheet1(:,1)-TripsDE10/1000).^2+(DE2011.data.Sheet1(:,1)-TripsDE11/1000).^2+(DE2012.data.Sheet1(:,1)-TripsDE12/1000).^2+(DE2013.data.Sheet1(:,1)-TripsDE13/1000).^2+(DE2014.data.Sheet1(:,1)-TripsDE14/1000).^2+(DE2015.data.Sheet1(:,1)-TripsDE15/1000).^2+(DE2016.data.Sheet1(:,1)-TripsDE16/1000).^2+(DE2017.data.Sheet1(:,1)-TripsDE17/1000).^2+(DE2018.data.Sheet1(:,1)-TripsDE18/1000).^2+(DE2019.data.Sheet1(:,1)-TripsDE19/1000).^2+(DE2020.data.Sheet1(:,1)-TripsDE20/1000).^2)/11;
MSPEM_DE=((DE2010.data.Sheet2(:,1)-TripsDE10/1000).^2+(DE2011.data.Sheet2(:,1)-TripsDE11/1000).^2+(DE2012.data.Sheet2(:,1)-TripsDE12/1000).^2+(DE2013.data.Sheet2(:,1)-TripsDE13/1000).^2+(DE2014.data.Sheet2(:,1)-TripsDE14/1000).^2+(DE2015.data.Sheet2(:,1)-TripsDE15/1000).^2+(DE2016.data.Sheet2(:,1)-TripsDE16/1000).^2+(DE2017.data.Sheet2(:,1)-TripsDE17/1000).^2+(DE2018.data.Sheet2(:,1)-TripsDE18/1000).^2+(DE2019.data.Sheet2(:,1)-TripsDE19/1000).^2+(DE2020.data.Sheet2(:,1)-TripsDE20/1000).^2)/11;
MSPEL_DE=((DE2010.data.Sheet3(:,1)-TripsDE10/1000).^2+(DE2011.data.Sheet3(:,1)-TripsDE11/1000).^2+(DE2012.data.Sheet3(:,1)-TripsDE12/1000).^2+(DE2013.data.Sheet3(:,1)-TripsDE13/1000).^2+(DE2014.data.Sheet3(:,1)-TripsDE14/1000).^2+(DE2015.data.Sheet3(:,1)-TripsDE15/1000).^2+(DE2016.data.Sheet3(:,1)-TripsDE16/1000).^2+(DE2017.data.Sheet3(:,1)-TripsDE17/1000).^2+(DE2018.data.Sheet3(:,1)-TripsDE18/1000).^2+(DE2019.data.Sheet3(:,1)-TripsDE19/1000).^2+(DE2020.data.Sheet3(:,1)-TripsDE20/1000).^2)/11;
MSPEC_DE=((DE2010.data.Sheet4(:,1)-TripsDE10/1000).^2+(DE2011.data.Sheet4(:,1)-TripsDE11/1000).^2+(DE2012.data.Sheet4(:,1)-TripsDE12/1000).^2+(DE2013.data.Sheet4(:,1)-TripsDE13/1000).^2+(DE2014.data.Sheet4(:,1)-TripsDE14/1000).^2+(DE2015.data.Sheet4(:,1)-TripsDE15/1000).^2+(DE2016.data.Sheet4(:,1)-TripsDE16/1000).^2+(DE2017.data.Sheet4(:,1)-TripsDE17/1000).^2+(DE2018.data.Sheet4(:,1)-TripsDE18/1000).^2+(DE2019.data.Sheet4(:,1)-TripsDE19/1000).^2+(DE2020.data.Sheet4(:,1)-TripsDE20/1000).^2)/11;

MSPEA_MD=((MD2010.data.Sheet1(:,1)-TripsMD10/1000).^2+(MD2011.data.Sheet1(:,1)-TripsMD11/1000).^2+(MD2012.data.Sheet1(:,1)-TripsMD12/1000).^2+(MD2013.data.Sheet1(:,1)-TripsMD13/1000).^2+(MD2014.data.Sheet1(:,1)-TripsMD14/1000).^2+(MD2015.data.Sheet1(:,1)-TripsMD15/1000).^2+(MD2016.data.Sheet1(:,1)-TripsMD16/1000).^2+(MD2017.data.Sheet1(:,1)-TripsMD17/1000).^2+(MD2018.data.Sheet1(:,1)-TripsMD18/1000).^2+(MD2019.data.Sheet1(:,1)-TripsMD19/1000).^2+(MD2020.data.Sheet1(:,1)-TripsMD20/1000).^2)/11;
MSPEM_MD=((MD2010.data.Sheet2(:,1)-TripsMD10/1000).^2+(MD2011.data.Sheet2(:,1)-TripsMD11/1000).^2+(MD2012.data.Sheet2(:,1)-TripsMD12/1000).^2+(MD2013.data.Sheet2(:,1)-TripsMD13/1000).^2+(MD2014.data.Sheet2(:,1)-TripsMD14/1000).^2+(MD2015.data.Sheet2(:,1)-TripsMD15/1000).^2+(MD2016.data.Sheet2(:,1)-TripsMD16/1000).^2+(MD2017.data.Sheet2(:,1)-TripsMD17/1000).^2+(MD2018.data.Sheet2(:,1)-TripsMD18/1000).^2+(MD2019.data.Sheet2(:,1)-TripsMD19/1000).^2+(MD2020.data.Sheet2(:,1)-TripsMD20/1000).^2)/11;
MSPEL_MD=((MD2010.data.Sheet3(:,1)-TripsMD10/1000).^2+(MD2011.data.Sheet3(:,1)-TripsMD11/1000).^2+(MD2012.data.Sheet3(:,1)-TripsMD12/1000).^2+(MD2013.data.Sheet3(:,1)-TripsMD13/1000).^2+(MD2014.data.Sheet3(:,1)-TripsMD14/1000).^2+(MD2015.data.Sheet3(:,1)-TripsMD15/1000).^2+(MD2016.data.Sheet3(:,1)-TripsMD16/1000).^2+(MD2017.data.Sheet3(:,1)-TripsMD17/1000).^2+(MD2018.data.Sheet3(:,1)-TripsMD18/1000).^2+(MD2019.data.Sheet3(:,1)-TripsMD19/1000).^2+(MD2020.data.Sheet3(:,1)-TripsMD20/1000).^2)/11;
MSPEC_MD=((MD2010.data.Sheet4(:,1)-TripsMD10/1000).^2+(MD2011.data.Sheet4(:,1)-TripsMD11/1000).^2+(MD2012.data.Sheet4(:,1)-TripsMD12/1000).^2+(MD2013.data.Sheet4(:,1)-TripsMD13/1000).^2+(MD2014.data.Sheet4(:,1)-TripsMD14/1000).^2+(MD2015.data.Sheet4(:,1)-TripsMD15/1000).^2+(MD2016.data.Sheet4(:,1)-TripsMD16/1000).^2+(MD2017.data.Sheet4(:,1)-TripsMD17/1000).^2+(MD2018.data.Sheet4(:,1)-TripsMD18/1000).^2+(MD2019.data.Sheet4(:,1)-TripsMD19/1000).^2+(MD2020.data.Sheet4(:,1)-TripsMD20/1000).^2)/11;

MSPEA_MA=((MA2010.data.Sheet1(:,1)-TripsMA10/1000).^2+(MA2011.data.Sheet1(:,1)-TripsMA11/1000).^2+(MA2012.data.Sheet1(:,1)-TripsMA12/1000).^2+(MA2013.data.Sheet1(:,1)-TripsMA13/1000).^2+(MA2014.data.Sheet1(:,1)-TripsMA14/1000).^2+(MA2015.data.Sheet1(:,1)-TripsMA15/1000).^2+(MA2016.data.Sheet1(:,1)-TripsMA16/1000).^2+(MA2017.data.Sheet1(:,1)-TripsMA17/1000).^2+(MA2018.data.Sheet1(:,1)-TripsMA18/1000).^2+(MA2019.data.Sheet1(:,1)-TripsMA19/1000).^2+(MA2020.data.Sheet1(:,1)-TripsMA20/1000).^2)/11;
MSPEM_MA=((MA2010.data.Sheet2(:,1)-TripsMA10/1000).^2+(MA2011.data.Sheet2(:,1)-TripsMA11/1000).^2+(MA2012.data.Sheet2(:,1)-TripsMA12/1000).^2+(MA2013.data.Sheet2(:,1)-TripsMA13/1000).^2+(MA2014.data.Sheet2(:,1)-TripsMA14/1000).^2+(MA2015.data.Sheet2(:,1)-TripsMA15/1000).^2+(MA2016.data.Sheet2(:,1)-TripsMA16/1000).^2+(MA2017.data.Sheet2(:,1)-TripsMA17/1000).^2+(MA2018.data.Sheet2(:,1)-TripsMA18/1000).^2+(MA2019.data.Sheet2(:,1)-TripsMA19/1000).^2+(MA2020.data.Sheet2(:,1)-TripsMA20/1000).^2)/11;
MSPEL_MA=((MA2010.data.Sheet3(:,1)-TripsMA10/1000).^2+(MA2011.data.Sheet3(:,1)-TripsMA11/1000).^2+(MA2012.data.Sheet3(:,1)-TripsMA12/1000).^2+(MA2013.data.Sheet3(:,1)-TripsMA13/1000).^2+(MA2014.data.Sheet3(:,1)-TripsMA14/1000).^2+(MA2015.data.Sheet3(:,1)-TripsMA15/1000).^2+(MA2016.data.Sheet3(:,1)-TripsMA16/1000).^2+(MA2017.data.Sheet3(:,1)-TripsMA17/1000).^2+(MA2018.data.Sheet3(:,1)-TripsMA18/1000).^2+(MA2019.data.Sheet3(:,1)-TripsMA19/1000).^2+(MA2020.data.Sheet3(:,1)-TripsMA20/1000).^2)/11;
MSPEC_MA=((MA2010.data.Sheet4(:,1)-TripsMA10/1000).^2+(MA2011.data.Sheet4(:,1)-TripsMA11/1000).^2+(MA2012.data.Sheet4(:,1)-TripsMA12/1000).^2+(MA2013.data.Sheet4(:,1)-TripsMA13/1000).^2+(MA2014.data.Sheet4(:,1)-TripsMA14/1000).^2+(MA2015.data.Sheet4(:,1)-TripsMA15/1000).^2+(MA2016.data.Sheet4(:,1)-TripsMA16/1000).^2+(MA2017.data.Sheet4(:,1)-TripsMA17/1000).^2+(MA2018.data.Sheet4(:,1)-TripsMA18/1000).^2+(MA2019.data.Sheet4(:,1)-TripsMA19/1000).^2+(MA2020.data.Sheet4(:,1)-TripsMA20/1000).^2)/11;

MSPEA_NJ=((NJ2010.data.Sheet1(:,1)-TripsNJ10/1000).^2+(NJ2011.data.Sheet1(:,1)-TripsNJ11/1000).^2+(NJ2012.data.Sheet1(:,1)-TripsNJ12/1000).^2+(NJ2013.data.Sheet1(:,1)-TripsNJ13/1000).^2+(NJ2014.data.Sheet1(:,1)-TripsNJ14/1000).^2+(NJ2015.data.Sheet1(:,1)-TripsNJ15/1000).^2+(NJ2016.data.Sheet1(:,1)-TripsNJ16/1000).^2+(NJ2017.data.Sheet1(:,1)-TripsNJ17/1000).^2+(NJ2018.data.Sheet1(:,1)-TripsNJ18/1000).^2+(NJ2019.data.Sheet1(:,1)-TripsNJ19/1000).^2+(NJ2020.data.Sheet1(:,1)-TripsNJ20/1000).^2)/11;
MSPEM_NJ=((NJ2010.data.Sheet2(:,1)-TripsNJ10/1000).^2+(NJ2011.data.Sheet2(:,1)-TripsNJ11/1000).^2+(NJ2012.data.Sheet2(:,1)-TripsNJ12/1000).^2+(NJ2013.data.Sheet2(:,1)-TripsNJ13/1000).^2+(NJ2014.data.Sheet2(:,1)-TripsNJ14/1000).^2+(NJ2015.data.Sheet2(:,1)-TripsNJ15/1000).^2+(NJ2016.data.Sheet2(:,1)-TripsNJ16/1000).^2+(NJ2017.data.Sheet2(:,1)-TripsNJ17/1000).^2+(NJ2018.data.Sheet2(:,1)-TripsNJ18/1000).^2+(NJ2019.data.Sheet2(:,1)-TripsNJ19/1000).^2+(NJ2020.data.Sheet2(:,1)-TripsNJ20/1000).^2)/11;
MSPEL_NJ=((NJ2010.data.Sheet3(:,1)-TripsNJ10/1000).^2+(NJ2011.data.Sheet3(:,1)-TripsNJ11/1000).^2+(NJ2012.data.Sheet3(:,1)-TripsNJ12/1000).^2+(NJ2013.data.Sheet3(:,1)-TripsNJ13/1000).^2+(NJ2014.data.Sheet3(:,1)-TripsNJ14/1000).^2+(NJ2015.data.Sheet3(:,1)-TripsNJ15/1000).^2+(NJ2016.data.Sheet3(:,1)-TripsNJ16/1000).^2+(NJ2017.data.Sheet3(:,1)-TripsNJ17/1000).^2+(NJ2018.data.Sheet3(:,1)-TripsNJ18/1000).^2+(NJ2019.data.Sheet3(:,1)-TripsNJ19/1000).^2+(NJ2020.data.Sheet3(:,1)-TripsNJ20/1000).^2)/11;
MSPEC_NJ=((NJ2010.data.Sheet4(:,1)-TripsNJ10/1000).^2+(NJ2011.data.Sheet4(:,1)-TripsNJ11/1000).^2+(NJ2012.data.Sheet4(:,1)-TripsNJ12/1000).^2+(NJ2013.data.Sheet4(:,1)-TripsNJ13/1000).^2+(NJ2014.data.Sheet4(:,1)-TripsNJ14/1000).^2+(NJ2015.data.Sheet4(:,1)-TripsNJ15/1000).^2+(NJ2016.data.Sheet4(:,1)-TripsNJ16/1000).^2+(NJ2017.data.Sheet4(:,1)-TripsNJ17/1000).^2+(NJ2018.data.Sheet4(:,1)-TripsNJ18/1000).^2+(NJ2019.data.Sheet4(:,1)-TripsNJ19/1000).^2+(NJ2020.data.Sheet4(:,1)-TripsNJ20/1000).^2)/11;

MSPEA_NY=((NY2010.data.Sheet1(:,1)-TripsNY10/1000).^2+(NY2011.data.Sheet1(:,1)-TripsNY11/1000).^2+(NY2012.data.Sheet1(:,1)-TripsNY12/1000).^2+(NY2013.data.Sheet1(:,1)-TripsNY13/1000).^2+(NY2014.data.Sheet1(:,1)-TripsNY14/1000).^2+(NY2015.data.Sheet1(:,1)-TripsNY15/1000).^2+(NY2016.data.Sheet1(:,1)-TripsNY16/1000).^2+(NY2017.data.Sheet1(:,1)-TripsNY17/1000).^2+(NY2018.data.Sheet1(:,1)-TripsNY18/1000).^2+(NY2019.data.Sheet1(:,1)-TripsNY19/1000).^2+(NY2020.data.Sheet1(:,1)-TripsNY20/1000).^2)/11;
MSPEM_NY=((NY2010.data.Sheet2(:,1)-TripsNY10/1000).^2+(NY2011.data.Sheet2(:,1)-TripsNY11/1000).^2+(NY2012.data.Sheet2(:,1)-TripsNY12/1000).^2+(NY2013.data.Sheet2(:,1)-TripsNY13/1000).^2+(NY2014.data.Sheet2(:,1)-TripsNY14/1000).^2+(NY2015.data.Sheet2(:,1)-TripsNY15/1000).^2+(NY2016.data.Sheet2(:,1)-TripsNY16/1000).^2+(NY2017.data.Sheet2(:,1)-TripsNY17/1000).^2+(NY2018.data.Sheet2(:,1)-TripsNY18/1000).^2+(NY2019.data.Sheet2(:,1)-TripsNY19/1000).^2+(NY2020.data.Sheet2(:,1)-TripsNY20/1000).^2)/11;
MSPEL_NY=((NY2010.data.Sheet3(:,1)-TripsNY10/1000).^2+(NY2011.data.Sheet3(:,1)-TripsNY11/1000).^2+(NY2012.data.Sheet3(:,1)-TripsNY12/1000).^2+(NY2013.data.Sheet3(:,1)-TripsNY13/1000).^2+(NY2014.data.Sheet3(:,1)-TripsNY14/1000).^2+(NY2015.data.Sheet3(:,1)-TripsNY15/1000).^2+(NY2016.data.Sheet3(:,1)-TripsNY16/1000).^2+(NY2017.data.Sheet3(:,1)-TripsNY17/1000).^2+(NY2018.data.Sheet3(:,1)-TripsNY18/1000).^2+(NY2019.data.Sheet3(:,1)-TripsNY19/1000).^2+(NY2020.data.Sheet3(:,1)-TripsNY20/1000).^2)/11;
MSPEC_NY=((NY2010.data.Sheet4(:,1)-TripsNY10/1000).^2+(NY2011.data.Sheet4(:,1)-TripsNY11/1000).^2+(NY2012.data.Sheet4(:,1)-TripsNY12/1000).^2+(NY2013.data.Sheet4(:,1)-TripsNY13/1000).^2+(NY2014.data.Sheet4(:,1)-TripsNY14/1000).^2+(NY2015.data.Sheet4(:,1)-TripsNY15/1000).^2+(NY2016.data.Sheet4(:,1)-TripsNY16/1000).^2+(NY2017.data.Sheet4(:,1)-TripsNY17/1000).^2+(NY2018.data.Sheet4(:,1)-TripsNY18/1000).^2+(NY2019.data.Sheet4(:,1)-TripsNY19/1000).^2+(NY2020.data.Sheet4(:,1)-TripsNY20/1000).^2)/11;

MSPEA_NC=((NC2010.data.Sheet1(:,1)-TripsNC10/1000).^2+(NC2011.data.Sheet1(:,1)-TripsNC11/1000).^2+(NC2012.data.Sheet1(:,1)-TripsNC12/1000).^2+(NC2013.data.Sheet1(:,1)-TripsNC13/1000).^2+(NC2014.data.Sheet1(:,1)-TripsNC14/1000).^2+(NC2015.data.Sheet1(:,1)-TripsNC15/1000).^2+(NC2016.data.Sheet1(:,1)-TripsNC16/1000).^2+(NC2017.data.Sheet1(:,1)-TripsNC17/1000).^2+(NC2018.data.Sheet1(:,1)-TripsNC18/1000).^2+(NC2019.data.Sheet1(:,1)-TripsNC19/1000).^2+(NC2020.data.Sheet1(:,1)-TripsNC20/1000).^2)/11;
MSPEM_NC=((NC2010.data.Sheet2(:,1)-TripsNC10/1000).^2+(NC2011.data.Sheet2(:,1)-TripsNC11/1000).^2+(NC2012.data.Sheet2(:,1)-TripsNC12/1000).^2+(NC2013.data.Sheet2(:,1)-TripsNC13/1000).^2+(NC2014.data.Sheet2(:,1)-TripsNC14/1000).^2+(NC2015.data.Sheet2(:,1)-TripsNC15/1000).^2+(NC2016.data.Sheet2(:,1)-TripsNC16/1000).^2+(NC2017.data.Sheet2(:,1)-TripsNC17/1000).^2+(NC2018.data.Sheet2(:,1)-TripsNC18/1000).^2+(NC2019.data.Sheet2(:,1)-TripsNC19/1000).^2+(NC2020.data.Sheet2(:,1)-TripsNC20/1000).^2)/11;
MSPEL_NC=((NC2010.data.Sheet3(:,1)-TripsNC10/1000).^2+(NC2011.data.Sheet3(:,1)-TripsNC11/1000).^2+(NC2012.data.Sheet3(:,1)-TripsNC12/1000).^2+(NC2013.data.Sheet3(:,1)-TripsNC13/1000).^2+(NC2014.data.Sheet3(:,1)-TripsNC14/1000).^2+(NC2015.data.Sheet3(:,1)-TripsNC15/1000).^2+(NC2016.data.Sheet3(:,1)-TripsNC16/1000).^2+(NC2017.data.Sheet3(:,1)-TripsNC17/1000).^2+(NC2018.data.Sheet3(:,1)-TripsNC18/1000).^2+(NC2019.data.Sheet3(:,1)-TripsNC19/1000).^2+(NC2020.data.Sheet3(:,1)-TripsNC20/1000).^2)/11;
MSPEC_NC=((NC2010.data.Sheet4(:,1)-TripsNC10/1000).^2+(NC2011.data.Sheet4(:,1)-TripsNC11/1000).^2+(NC2012.data.Sheet4(:,1)-TripsNC12/1000).^2+(NC2013.data.Sheet4(:,1)-TripsNC13/1000).^2+(NC2014.data.Sheet4(:,1)-TripsNC14/1000).^2+(NC2015.data.Sheet4(:,1)-TripsNC15/1000).^2+(NC2016.data.Sheet4(:,1)-TripsNC16/1000).^2+(NC2017.data.Sheet4(:,1)-TripsNC17/1000).^2+(NC2018.data.Sheet4(:,1)-TripsNC18/1000).^2+(NC2019.data.Sheet4(:,1)-TripsNC19/1000).^2+(NC2020.data.Sheet4(:,1)-TripsNC20/1000).^2)/11;

MSPEA_RI=((RI2010.data.Sheet1(:,1)-TripsRI10/1000).^2+(RI2011.data.Sheet1(:,1)-TripsRI11/1000).^2+(RI2012.data.Sheet1(:,1)-TripsRI12/1000).^2+(RI2013.data.Sheet1(:,1)-TripsRI13/1000).^2+(RI2014.data.Sheet1(:,1)-TripsRI14/1000).^2+(RI2015.data.Sheet1(:,1)-TripsRI15/1000).^2+(RI2016.data.Sheet1(:,1)-TripsRI16/1000).^2+(RI2017.data.Sheet1(:,1)-TripsRI17/1000).^2+(RI2018.data.Sheet1(:,1)-TripsRI18/1000).^2+(RI2019.data.Sheet1(:,1)-TripsRI19/1000).^2+(RI2020.data.Sheet1(:,1)-TripsRI20/1000).^2)/11;
MSPEM_RI=((RI2010.data.Sheet2(:,1)-TripsRI10/1000).^2+(RI2011.data.Sheet2(:,1)-TripsRI11/1000).^2+(RI2012.data.Sheet2(:,1)-TripsRI12/1000).^2+(RI2013.data.Sheet2(:,1)-TripsRI13/1000).^2+(RI2014.data.Sheet2(:,1)-TripsRI14/1000).^2+(RI2015.data.Sheet2(:,1)-TripsRI15/1000).^2+(RI2016.data.Sheet2(:,1)-TripsRI16/1000).^2+(RI2017.data.Sheet2(:,1)-TripsRI17/1000).^2+(RI2018.data.Sheet2(:,1)-TripsRI18/1000).^2+(RI2019.data.Sheet2(:,1)-TripsRI19/1000).^2+(RI2020.data.Sheet2(:,1)-TripsRI20/1000).^2)/11;
MSPEL_RI=((RI2010.data.Sheet3(:,1)-TripsRI10/1000).^2+(RI2011.data.Sheet3(:,1)-TripsRI11/1000).^2+(RI2012.data.Sheet3(:,1)-TripsRI12/1000).^2+(RI2013.data.Sheet3(:,1)-TripsRI13/1000).^2+(RI2014.data.Sheet3(:,1)-TripsRI14/1000).^2+(RI2015.data.Sheet3(:,1)-TripsRI15/1000).^2+(RI2016.data.Sheet3(:,1)-TripsRI16/1000).^2+(RI2017.data.Sheet3(:,1)-TripsRI17/1000).^2+(RI2018.data.Sheet3(:,1)-TripsRI18/1000).^2+(RI2019.data.Sheet3(:,1)-TripsRI19/1000).^2+(RI2020.data.Sheet3(:,1)-TripsRI20/1000).^2)/11;
MSPEC_RI=((RI2010.data.Sheet4(:,1)-TripsRI10/1000).^2+(RI2011.data.Sheet4(:,1)-TripsRI11/1000).^2+(RI2012.data.Sheet4(:,1)-TripsRI12/1000).^2+(RI2013.data.Sheet4(:,1)-TripsRI13/1000).^2+(RI2014.data.Sheet4(:,1)-TripsRI14/1000).^2+(RI2015.data.Sheet4(:,1)-TripsRI15/1000).^2+(RI2016.data.Sheet4(:,1)-TripsRI16/1000).^2+(RI2017.data.Sheet4(:,1)-TripsRI17/1000).^2+(RI2018.data.Sheet4(:,1)-TripsRI18/1000).^2+(RI2019.data.Sheet4(:,1)-TripsRI19/1000).^2+(RI2020.data.Sheet4(:,1)-TripsRI20/1000).^2)/11;

MSPEA_VA=((VA2010.data.Sheet1(:,1)-TripsVA10/1000).^2+(VA2011.data.Sheet1(:,1)-TripsVA11/1000).^2+(VA2012.data.Sheet1(:,1)-TripsVA12/1000).^2+(VA2013.data.Sheet1(:,1)-TripsVA13/1000).^2+(VA2014.data.Sheet1(:,1)-TripsVA14/1000).^2+(VA2015.data.Sheet1(:,1)-TripsVA15/1000).^2+(VA2016.data.Sheet1(:,1)-TripsVA16/1000).^2+(VA2017.data.Sheet1(:,1)-TripsVA17/1000).^2+(VA2018.data.Sheet1(:,1)-TripsVA18/1000).^2+(VA2019.data.Sheet1(:,1)-TripsVA19/1000).^2+(VA2020.data.Sheet1(:,1)-TripsVA20/1000).^2)/11;
MSPEM_VA=((VA2010.data.Sheet2(:,1)-TripsVA10/1000).^2+(VA2011.data.Sheet2(:,1)-TripsVA11/1000).^2+(VA2012.data.Sheet2(:,1)-TripsVA12/1000).^2+(VA2013.data.Sheet2(:,1)-TripsVA13/1000).^2+(VA2014.data.Sheet2(:,1)-TripsVA14/1000).^2+(VA2015.data.Sheet2(:,1)-TripsVA15/1000).^2+(VA2016.data.Sheet2(:,1)-TripsVA16/1000).^2+(VA2017.data.Sheet2(:,1)-TripsVA17/1000).^2+(VA2018.data.Sheet2(:,1)-TripsVA18/1000).^2+(VA2019.data.Sheet2(:,1)-TripsVA19/1000).^2+(VA2020.data.Sheet2(:,1)-TripsVA20/1000).^2)/11;
MSPEL_VA=((VA2010.data.Sheet3(:,1)-TripsVA10/1000).^2+(VA2011.data.Sheet3(:,1)-TripsVA11/1000).^2+(VA2012.data.Sheet3(:,1)-TripsVA12/1000).^2+(VA2013.data.Sheet3(:,1)-TripsVA13/1000).^2+(VA2014.data.Sheet3(:,1)-TripsVA14/1000).^2+(VA2015.data.Sheet3(:,1)-TripsVA15/1000).^2+(VA2016.data.Sheet3(:,1)-TripsVA16/1000).^2+(VA2017.data.Sheet3(:,1)-TripsVA17/1000).^2+(VA2018.data.Sheet3(:,1)-TripsVA18/1000).^2+(VA2019.data.Sheet3(:,1)-TripsVA19/1000).^2+(VA2020.data.Sheet3(:,1)-TripsVA20/1000).^2)/11;
MSPEC_VA=((VA2010.data.Sheet4(:,1)-TripsVA10/1000).^2+(VA2011.data.Sheet4(:,1)-TripsVA11/1000).^2+(VA2012.data.Sheet4(:,1)-TripsVA12/1000).^2+(VA2013.data.Sheet4(:,1)-TripsVA13/1000).^2+(VA2014.data.Sheet4(:,1)-TripsVA14/1000).^2+(VA2015.data.Sheet4(:,1)-TripsVA15/1000).^2+(VA2016.data.Sheet4(:,1)-TripsVA16/1000).^2+(VA2017.data.Sheet4(:,1)-TripsVA17/1000).^2+(VA2018.data.Sheet4(:,1)-TripsVA18/1000).^2+(VA2019.data.Sheet4(:,1)-TripsVA19/1000).^2+(VA2020.data.Sheet4(:,1)-TripsVA20/1000).^2)/11;

% Landings
MSPEFishA_CT=((CT2010.data.Sheet1(:,4)-FishCT10/1000).^2+(CT2011.data.Sheet1(:,4)-FishCT11/1000).^2+(CT2012.data.Sheet1(:,4)-FishCT12/1000).^2+(CT2013.data.Sheet1(:,4)-FishCT13/1000).^2+(CT2014.data.Sheet1(:,4)-FishCT14/1000).^2+(CT2015.data.Sheet1(:,4)-FishCT15/1000).^2+(CT2016.data.Sheet1(:,4)-FishCT16/1000).^2+(CT2017.data.Sheet1(:,4)-FishCT17/1000).^2+(CT2018.data.Sheet1(:,4)-FishCT18/1000).^2+(CT2019.data.Sheet1(:,4)-FishCT19/1000).^2+(CT2020.data.Sheet1(:,4)-FishCT20/1000).^2)/11;
MSPEFishM_CT=((CT2010.data.Sheet2(:,4)-FishCT10/1000).^2+(CT2011.data.Sheet2(:,4)-FishCT11/1000).^2+(CT2012.data.Sheet2(:,4)-FishCT12/1000).^2+(CT2013.data.Sheet2(:,4)-FishCT13/1000).^2+(CT2014.data.Sheet2(:,4)-FishCT14/1000).^2+(CT2015.data.Sheet2(:,4)-FishCT15/1000).^2+(CT2016.data.Sheet2(:,4)-FishCT16/1000).^2+(CT2017.data.Sheet2(:,4)-FishCT17/1000).^2+(CT2018.data.Sheet2(:,4)-FishCT18/1000).^2+(CT2019.data.Sheet2(:,4)-FishCT19/1000).^2+(CT2020.data.Sheet2(:,4)-FishCT20/1000).^2)/11;
MSPEFishL_CT=((CT2010.data.Sheet3(:,4)-FishCT10/1000).^2+(CT2011.data.Sheet3(:,4)-FishCT11/1000).^2+(CT2012.data.Sheet3(:,4)-FishCT12/1000).^2+(CT2013.data.Sheet3(:,4)-FishCT13/1000).^2+(CT2014.data.Sheet3(:,4)-FishCT14/1000).^2+(CT2015.data.Sheet3(:,4)-FishCT15/1000).^2+(CT2016.data.Sheet3(:,4)-FishCT16/1000).^2+(CT2017.data.Sheet3(:,4)-FishCT17/1000).^2+(CT2018.data.Sheet3(:,4)-FishCT18/1000).^2+(CT2019.data.Sheet3(:,4)-FishCT19/1000).^2+(CT2020.data.Sheet3(:,4)-FishCT20/1000).^2)/11;
MSPEFishC_CT=((CT2010.data.Sheet4(:,4)-FishCT10/1000).^2+(CT2011.data.Sheet4(:,4)-FishCT11/1000).^2+(CT2012.data.Sheet4(:,4)-FishCT12/1000).^2+(CT2013.data.Sheet4(:,4)-FishCT13/1000).^2+(CT2014.data.Sheet4(:,4)-FishCT14/1000).^2+(CT2015.data.Sheet4(:,4)-FishCT15/1000).^2+(CT2016.data.Sheet4(:,4)-FishCT16/1000).^2+(CT2017.data.Sheet4(:,4)-FishCT17/1000).^2+(CT2018.data.Sheet4(:,4)-FishCT18/1000).^2+(CT2019.data.Sheet4(:,4)-FishCT19/1000).^2+(CT2020.data.Sheet4(:,4)-FishCT20/1000).^2)/11;

MSPEFishA_DE=((DE2010.data.Sheet1(:,4)-FishDE10/1000).^2+(DE2011.data.Sheet1(:,4)-FishDE11/1000).^2+(DE2012.data.Sheet1(:,4)-FishDE12/1000).^2+(DE2013.data.Sheet1(:,4)-FishDE13/1000).^2+(DE2014.data.Sheet1(:,4)-FishDE14/1000).^2+(DE2015.data.Sheet1(:,4)-FishDE15/1000).^2+(DE2016.data.Sheet1(:,4)-FishDE16/1000).^2+(DE2017.data.Sheet1(:,4)-FishDE17/1000).^2+(DE2018.data.Sheet1(:,4)-FishDE18/1000).^2+(DE2019.data.Sheet1(:,4)-FishDE19/1000).^2+(DE2020.data.Sheet1(:,4)-FishDE20/1000).^2)/11;
MSPEFishM_DE=((DE2010.data.Sheet2(:,4)-FishDE10/1000).^2+(DE2011.data.Sheet2(:,4)-FishDE11/1000).^2+(DE2012.data.Sheet2(:,4)-FishDE12/1000).^2+(DE2013.data.Sheet2(:,4)-FishDE13/1000).^2+(DE2014.data.Sheet2(:,4)-FishDE14/1000).^2+(DE2015.data.Sheet2(:,4)-FishDE15/1000).^2+(DE2016.data.Sheet2(:,4)-FishDE16/1000).^2+(DE2017.data.Sheet2(:,4)-FishDE17/1000).^2+(DE2018.data.Sheet2(:,4)-FishDE18/1000).^2+(DE2019.data.Sheet2(:,4)-FishDE19/1000).^2+(DE2020.data.Sheet2(:,4)-FishDE20/1000).^2)/11;
MSPEFishL_DE=((DE2010.data.Sheet3(:,4)-FishDE10/1000).^2+(DE2011.data.Sheet3(:,4)-FishDE11/1000).^2+(DE2012.data.Sheet3(:,4)-FishDE12/1000).^2+(DE2013.data.Sheet3(:,4)-FishDE13/1000).^2+(DE2014.data.Sheet3(:,4)-FishDE14/1000).^2+(DE2015.data.Sheet3(:,4)-FishDE15/1000).^2+(DE2016.data.Sheet3(:,4)-FishDE16/1000).^2+(DE2017.data.Sheet3(:,4)-FishDE17/1000).^2+(DE2018.data.Sheet3(:,4)-FishDE18/1000).^2+(DE2019.data.Sheet3(:,4)-FishDE19/1000).^2+(DE2020.data.Sheet3(:,4)-FishDE20/1000).^2)/11;
MSPEFishC_DE=((DE2010.data.Sheet4(:,4)-FishDE10/1000).^2+(DE2011.data.Sheet4(:,4)-FishDE11/1000).^2+(DE2012.data.Sheet4(:,4)-FishDE12/1000).^2+(DE2013.data.Sheet4(:,4)-FishDE13/1000).^2+(DE2014.data.Sheet4(:,4)-FishDE14/1000).^2+(DE2015.data.Sheet4(:,4)-FishDE15/1000).^2+(DE2016.data.Sheet4(:,4)-FishDE16/1000).^2+(DE2017.data.Sheet4(:,4)-FishDE17/1000).^2+(DE2018.data.Sheet4(:,4)-FishDE18/1000).^2+(DE2019.data.Sheet4(:,4)-FishDE19/1000).^2+(DE2020.data.Sheet4(:,4)-FishDE20/1000).^2)/11;

MSPEFishA_MD=((MD2010.data.Sheet1(:,4)-FishMD10/1000).^2+(MD2011.data.Sheet1(:,4)-FishMD11/1000).^2+(MD2012.data.Sheet1(:,4)-FishMD12/1000).^2+(MD2013.data.Sheet1(:,4)-FishMD13/1000).^2+(MD2014.data.Sheet1(:,4)-FishMD14/1000).^2+(MD2015.data.Sheet1(:,4)-FishMD15/1000).^2+(MD2016.data.Sheet1(:,4)-FishMD16/1000).^2+(MD2017.data.Sheet1(:,4)-FishMD17/1000).^2+(MD2018.data.Sheet1(:,4)-FishMD18/1000).^2+(MD2019.data.Sheet1(:,4)-FishMD19/1000).^2+(MD2020.data.Sheet1(:,4)-FishMD20/1000).^2)/11;
MSPEFishM_MD=((MD2010.data.Sheet2(:,4)-FishMD10/1000).^2+(MD2011.data.Sheet2(:,4)-FishMD11/1000).^2+(MD2012.data.Sheet2(:,4)-FishMD12/1000).^2+(MD2013.data.Sheet2(:,4)-FishMD13/1000).^2+(MD2014.data.Sheet2(:,4)-FishMD14/1000).^2+(MD2015.data.Sheet2(:,4)-FishMD15/1000).^2+(MD2016.data.Sheet2(:,4)-FishMD16/1000).^2+(MD2017.data.Sheet2(:,4)-FishMD17/1000).^2+(MD2018.data.Sheet2(:,4)-FishMD18/1000).^2+(MD2019.data.Sheet2(:,4)-FishMD19/1000).^2+(MD2020.data.Sheet2(:,4)-FishMD20/1000).^2)/11;
MSPEFishL_MD=((MD2010.data.Sheet3(:,4)-FishMD10/1000).^2+(MD2011.data.Sheet3(:,4)-FishMD11/1000).^2+(MD2012.data.Sheet3(:,4)-FishMD12/1000).^2+(MD2013.data.Sheet3(:,4)-FishMD13/1000).^2+(MD2014.data.Sheet3(:,4)-FishMD14/1000).^2+(MD2015.data.Sheet3(:,4)-FishMD15/1000).^2+(MD2016.data.Sheet3(:,4)-FishMD16/1000).^2+(MD2017.data.Sheet3(:,4)-FishMD17/1000).^2+(MD2018.data.Sheet3(:,4)-FishMD18/1000).^2+(MD2019.data.Sheet3(:,4)-FishMD19/1000).^2+(MD2020.data.Sheet3(:,4)-FishMD20/1000).^2)/11;
MSPEFishC_MD=((MD2010.data.Sheet4(:,4)-FishMD10/1000).^2+(MD2011.data.Sheet4(:,4)-FishMD11/1000).^2+(MD2012.data.Sheet4(:,4)-FishMD12/1000).^2+(MD2013.data.Sheet4(:,4)-FishMD13/1000).^2+(MD2014.data.Sheet4(:,4)-FishMD14/1000).^2+(MD2015.data.Sheet4(:,4)-FishMD15/1000).^2+(MD2016.data.Sheet4(:,4)-FishMD16/1000).^2+(MD2017.data.Sheet4(:,4)-FishMD17/1000).^2+(MD2018.data.Sheet4(:,4)-FishMD18/1000).^2+(MD2019.data.Sheet4(:,4)-FishMD19/1000).^2+(MD2020.data.Sheet4(:,4)-FishMD20/1000).^2)/11;

MSPEFishA_MA=((MA2010.data.Sheet1(:,4)-FishMA10/1000).^2+(MA2011.data.Sheet1(:,4)-FishMA11/1000).^2+(MA2012.data.Sheet1(:,4)-FishMA12/1000).^2+(MA2013.data.Sheet1(:,4)-FishMA13/1000).^2+(MA2014.data.Sheet1(:,4)-FishMA14/1000).^2+(MA2015.data.Sheet1(:,4)-FishMA15/1000).^2+(MA2016.data.Sheet1(:,4)-FishMA16/1000).^2+(MA2017.data.Sheet1(:,4)-FishMA17/1000).^2+(MA2018.data.Sheet1(:,4)-FishMA18/1000).^2+(MA2019.data.Sheet1(:,4)-FishMA19/1000).^2+(MA2020.data.Sheet1(:,4)-FishMA20/1000).^2)/11;
MSPEFishM_MA=((MA2010.data.Sheet2(:,4)-FishMA10/1000).^2+(MA2011.data.Sheet2(:,4)-FishMA11/1000).^2+(MA2012.data.Sheet2(:,4)-FishMA12/1000).^2+(MA2013.data.Sheet2(:,4)-FishMA13/1000).^2+(MA2014.data.Sheet2(:,4)-FishMA14/1000).^2+(MA2015.data.Sheet2(:,4)-FishMA15/1000).^2+(MA2016.data.Sheet2(:,4)-FishMA16/1000).^2+(MA2017.data.Sheet2(:,4)-FishMA17/1000).^2+(MA2018.data.Sheet2(:,4)-FishMA18/1000).^2+(MA2019.data.Sheet2(:,4)-FishMA19/1000).^2+(MA2020.data.Sheet2(:,4)-FishMA20/1000).^2)/11;
MSPEFishL_MA=((MA2010.data.Sheet3(:,4)-FishMA10/1000).^2+(MA2011.data.Sheet3(:,4)-FishMA11/1000).^2+(MA2012.data.Sheet3(:,4)-FishMA12/1000).^2+(MA2013.data.Sheet3(:,4)-FishMA13/1000).^2+(MA2014.data.Sheet3(:,4)-FishMA14/1000).^2+(MA2015.data.Sheet3(:,4)-FishMA15/1000).^2+(MA2016.data.Sheet3(:,4)-FishMA16/1000).^2+(MA2017.data.Sheet3(:,4)-FishMA17/1000).^2+(MA2018.data.Sheet3(:,4)-FishMA18/1000).^2+(MA2019.data.Sheet3(:,4)-FishMA19/1000).^2+(MA2020.data.Sheet3(:,4)-FishMA20/1000).^2)/11;
MSPEFishC_MA=((MA2010.data.Sheet4(:,4)-FishMA10/1000).^2+(MA2011.data.Sheet4(:,4)-FishMA11/1000).^2+(MA2012.data.Sheet4(:,4)-FishMA12/1000).^2+(MA2013.data.Sheet4(:,4)-FishMA13/1000).^2+(MA2014.data.Sheet4(:,4)-FishMA14/1000).^2+(MA2015.data.Sheet4(:,4)-FishMA15/1000).^2+(MA2016.data.Sheet4(:,4)-FishMA16/1000).^2+(MA2017.data.Sheet4(:,4)-FishMA17/1000).^2+(MA2018.data.Sheet4(:,4)-FishMA18/1000).^2+(MA2019.data.Sheet4(:,4)-FishMA19/1000).^2+(MA2020.data.Sheet4(:,4)-FishMA20/1000).^2)/11;

MSPEFishA_NJ=((NJ2010.data.Sheet1(:,4)-FishNJ10/1000).^2+(NJ2011.data.Sheet1(:,4)-FishNJ11/1000).^2+(NJ2012.data.Sheet1(:,4)-FishNJ12/1000).^2+(NJ2013.data.Sheet1(:,4)-FishNJ13/1000).^2+(NJ2014.data.Sheet1(:,4)-FishNJ14/1000).^2+(NJ2015.data.Sheet1(:,4)-FishNJ15/1000).^2+(NJ2016.data.Sheet1(:,4)-FishNJ16/1000).^2+(NJ2017.data.Sheet1(:,4)-FishNJ17/1000).^2+(NJ2018.data.Sheet1(:,4)-FishNJ18/1000).^2+(NJ2019.data.Sheet1(:,4)-FishNJ19/1000).^2+(NJ2020.data.Sheet1(:,4)-FishNJ20/1000).^2)/11;
MSPEFishM_NJ=((NJ2010.data.Sheet2(:,4)-FishNJ10/1000).^2+(NJ2011.data.Sheet2(:,4)-FishNJ11/1000).^2+(NJ2012.data.Sheet2(:,4)-FishNJ12/1000).^2+(NJ2013.data.Sheet2(:,4)-FishNJ13/1000).^2+(NJ2014.data.Sheet2(:,4)-FishNJ14/1000).^2+(NJ2015.data.Sheet2(:,4)-FishNJ15/1000).^2+(NJ2016.data.Sheet2(:,4)-FishNJ16/1000).^2+(NJ2017.data.Sheet2(:,4)-FishNJ17/1000).^2+(NJ2018.data.Sheet2(:,4)-FishNJ18/1000).^2+(NJ2019.data.Sheet2(:,4)-FishNJ19/1000).^2+(NJ2020.data.Sheet2(:,4)-FishNJ20/1000).^2)/11;
MSPEFishL_NJ=((NJ2010.data.Sheet3(:,4)-FishNJ10/1000).^2+(NJ2011.data.Sheet3(:,4)-FishNJ11/1000).^2+(NJ2012.data.Sheet3(:,4)-FishNJ12/1000).^2+(NJ2013.data.Sheet3(:,4)-FishNJ13/1000).^2+(NJ2014.data.Sheet3(:,4)-FishNJ14/1000).^2+(NJ2015.data.Sheet3(:,4)-FishNJ15/1000).^2+(NJ2016.data.Sheet3(:,4)-FishNJ16/1000).^2+(NJ2017.data.Sheet3(:,4)-FishNJ17/1000).^2+(NJ2018.data.Sheet3(:,4)-FishNJ18/1000).^2+(NJ2019.data.Sheet3(:,4)-FishNJ19/1000).^2+(NJ2020.data.Sheet3(:,4)-FishNJ20/1000).^2)/11;
MSPEFishC_NJ=((NJ2010.data.Sheet4(:,4)-FishNJ10/1000).^2+(NJ2011.data.Sheet4(:,4)-FishNJ11/1000).^2+(NJ2012.data.Sheet4(:,4)-FishNJ12/1000).^2+(NJ2013.data.Sheet4(:,4)-FishNJ13/1000).^2+(NJ2014.data.Sheet4(:,4)-FishNJ14/1000).^2+(NJ2015.data.Sheet4(:,4)-FishNJ15/1000).^2+(NJ2016.data.Sheet4(:,4)-FishNJ16/1000).^2+(NJ2017.data.Sheet4(:,4)-FishNJ17/1000).^2+(NJ2018.data.Sheet4(:,4)-FishNJ18/1000).^2+(NJ2019.data.Sheet4(:,4)-FishNJ19/1000).^2+(NJ2020.data.Sheet4(:,4)-FishNJ20/1000).^2)/11;

MSPEFishA_NY=((NY2010.data.Sheet1(:,4)-FishNY10/1000).^2+(NY2011.data.Sheet1(:,4)-FishNY11/1000).^2+(NY2012.data.Sheet1(:,4)-FishNY12/1000).^2+(NY2013.data.Sheet1(:,4)-FishNY13/1000).^2+(NY2014.data.Sheet1(:,4)-FishNY14/1000).^2+(NY2015.data.Sheet1(:,4)-FishNY15/1000).^2+(NY2016.data.Sheet1(:,4)-FishNY16/1000).^2+(NY2017.data.Sheet1(:,4)-FishNY17/1000).^2+(NY2018.data.Sheet1(:,4)-FishNY18/1000).^2+(NY2019.data.Sheet1(:,4)-FishNY19/1000).^2+(NY2020.data.Sheet1(:,4)-FishNY20/1000).^2)/11;
MSPEFishM_NY=((NY2010.data.Sheet2(:,4)-FishNY10/1000).^2+(NY2011.data.Sheet2(:,4)-FishNY11/1000).^2+(NY2012.data.Sheet2(:,4)-FishNY12/1000).^2+(NY2013.data.Sheet2(:,4)-FishNY13/1000).^2+(NY2014.data.Sheet2(:,4)-FishNY14/1000).^2+(NY2015.data.Sheet2(:,4)-FishNY15/1000).^2+(NY2016.data.Sheet2(:,4)-FishNY16/1000).^2+(NY2017.data.Sheet2(:,4)-FishNY17/1000).^2+(NY2018.data.Sheet2(:,4)-FishNY18/1000).^2+(NY2019.data.Sheet2(:,4)-FishNY19/1000).^2+(NY2020.data.Sheet2(:,4)-FishNY20/1000).^2)/11;
MSPEFishL_NY=((NY2010.data.Sheet3(:,4)-FishNY10/1000).^2+(NY2011.data.Sheet3(:,4)-FishNY11/1000).^2+(NY2012.data.Sheet3(:,4)-FishNY12/1000).^2+(NY2013.data.Sheet3(:,4)-FishNY13/1000).^2+(NY2014.data.Sheet3(:,4)-FishNY14/1000).^2+(NY2015.data.Sheet3(:,4)-FishNY15/1000).^2+(NY2016.data.Sheet3(:,4)-FishNY16/1000).^2+(NY2017.data.Sheet3(:,4)-FishNY17/1000).^2+(NY2018.data.Sheet3(:,4)-FishNY18/1000).^2+(NY2019.data.Sheet3(:,4)-FishNY19/1000).^2+(NY2020.data.Sheet3(:,4)-FishNY20/1000).^2)/11;
MSPEFishC_NY=((NY2010.data.Sheet4(:,4)-FishNY10/1000).^2+(NY2011.data.Sheet4(:,4)-FishNY11/1000).^2+(NY2012.data.Sheet4(:,4)-FishNY12/1000).^2+(NY2013.data.Sheet4(:,4)-FishNY13/1000).^2+(NY2014.data.Sheet4(:,4)-FishNY14/1000).^2+(NY2015.data.Sheet4(:,4)-FishNY15/1000).^2+(NY2016.data.Sheet4(:,4)-FishNY16/1000).^2+(NY2017.data.Sheet4(:,4)-FishNY17/1000).^2+(NY2018.data.Sheet4(:,4)-FishNY18/1000).^2+(NY2019.data.Sheet4(:,4)-FishNY19/1000).^2+(NY2020.data.Sheet4(:,4)-FishNY20/1000).^2)/11;

MSPEFishA_NC=((NC2010.data.Sheet1(:,4)-FishNC10/1000).^2+(NC2011.data.Sheet1(:,4)-FishNC11/1000).^2+(NC2012.data.Sheet1(:,4)-FishNC12/1000).^2+(NC2013.data.Sheet1(:,4)-FishNC13/1000).^2+(NC2014.data.Sheet1(:,4)-FishNC14/1000).^2+(NC2015.data.Sheet1(:,4)-FishNC15/1000).^2+(NC2016.data.Sheet1(:,4)-FishNC16/1000).^2+(NC2017.data.Sheet1(:,4)-FishNC17/1000).^2+(NC2018.data.Sheet1(:,4)-FishNC18/1000).^2+(NC2019.data.Sheet1(:,4)-FishNC19/1000).^2+(NC2020.data.Sheet1(:,4)-FishNC20/1000).^2)/11;
MSPEFishM_NC=((NC2010.data.Sheet2(:,4)-FishNC10/1000).^2+(NC2011.data.Sheet2(:,4)-FishNC11/1000).^2+(NC2012.data.Sheet2(:,4)-FishNC12/1000).^2+(NC2013.data.Sheet2(:,4)-FishNC13/1000).^2+(NC2014.data.Sheet2(:,4)-FishNC14/1000).^2+(NC2015.data.Sheet2(:,4)-FishNC15/1000).^2+(NC2016.data.Sheet2(:,4)-FishNC16/1000).^2+(NC2017.data.Sheet2(:,4)-FishNC17/1000).^2+(NC2018.data.Sheet2(:,4)-FishNC18/1000).^2+(NC2019.data.Sheet2(:,4)-FishNC19/1000).^2+(NC2020.data.Sheet2(:,4)-FishNC20/1000).^2)/11;
MSPEFishL_NC=((NC2010.data.Sheet3(:,4)-FishNC10/1000).^2+(NC2011.data.Sheet3(:,4)-FishNC11/1000).^2+(NC2012.data.Sheet3(:,4)-FishNC12/1000).^2+(NC2013.data.Sheet3(:,4)-FishNC13/1000).^2+(NC2014.data.Sheet3(:,4)-FishNC14/1000).^2+(NC2015.data.Sheet3(:,4)-FishNC15/1000).^2+(NC2016.data.Sheet3(:,4)-FishNC16/1000).^2+(NC2017.data.Sheet3(:,4)-FishNC17/1000).^2+(NC2018.data.Sheet3(:,4)-FishNC18/1000).^2+(NC2019.data.Sheet3(:,4)-FishNC19/1000).^2+(NC2020.data.Sheet3(:,4)-FishNC20/1000).^2)/11;
MSPEFishC_NC=((NC2010.data.Sheet4(:,4)-FishNC10/1000).^2+(NC2011.data.Sheet4(:,4)-FishNC11/1000).^2+(NC2012.data.Sheet4(:,4)-FishNC12/1000).^2+(NC2013.data.Sheet4(:,4)-FishNC13/1000).^2+(NC2014.data.Sheet4(:,4)-FishNC14/1000).^2+(NC2015.data.Sheet4(:,4)-FishNC15/1000).^2+(NC2016.data.Sheet4(:,4)-FishNC16/1000).^2+(NC2017.data.Sheet4(:,4)-FishNC17/1000).^2+(NC2018.data.Sheet4(:,4)-FishNC18/1000).^2+(NC2019.data.Sheet4(:,4)-FishNC19/1000).^2+(NC2020.data.Sheet4(:,4)-FishNC20/1000).^2)/11;

MSPEFishA_RI=((RI2010.data.Sheet1(:,4)-FishRI10/1000).^2+(RI2011.data.Sheet1(:,4)-FishRI11/1000).^2+(RI2012.data.Sheet1(:,4)-FishRI12/1000).^2+(RI2013.data.Sheet1(:,4)-FishRI13/1000).^2+(RI2014.data.Sheet1(:,4)-FishRI14/1000).^2+(RI2015.data.Sheet1(:,4)-FishRI15/1000).^2+(RI2016.data.Sheet1(:,4)-FishRI16/1000).^2+(RI2017.data.Sheet1(:,4)-FishRI17/1000).^2+(RI2018.data.Sheet1(:,4)-FishRI18/1000).^2+(RI2019.data.Sheet1(:,4)-FishRI19/1000).^2+(RI2020.data.Sheet1(:,4)-FishRI20/1000).^2)/11;
MSPEFishM_RI=((RI2010.data.Sheet2(:,4)-FishRI10/1000).^2+(RI2011.data.Sheet2(:,4)-FishRI11/1000).^2+(RI2012.data.Sheet2(:,4)-FishRI12/1000).^2+(RI2013.data.Sheet2(:,4)-FishRI13/1000).^2+(RI2014.data.Sheet2(:,4)-FishRI14/1000).^2+(RI2015.data.Sheet2(:,4)-FishRI15/1000).^2+(RI2016.data.Sheet2(:,4)-FishRI16/1000).^2+(RI2017.data.Sheet2(:,4)-FishRI17/1000).^2+(RI2018.data.Sheet2(:,4)-FishRI18/1000).^2+(RI2019.data.Sheet2(:,4)-FishRI19/1000).^2+(RI2020.data.Sheet2(:,4)-FishRI20/1000).^2)/11;
MSPEFishL_RI=((RI2010.data.Sheet3(:,4)-FishRI10/1000).^2+(RI2011.data.Sheet3(:,4)-FishRI11/1000).^2+(RI2012.data.Sheet3(:,4)-FishRI12/1000).^2+(RI2013.data.Sheet3(:,4)-FishRI13/1000).^2+(RI2014.data.Sheet3(:,4)-FishRI14/1000).^2+(RI2015.data.Sheet3(:,4)-FishRI15/1000).^2+(RI2016.data.Sheet3(:,4)-FishRI16/1000).^2+(RI2017.data.Sheet3(:,4)-FishRI17/1000).^2+(RI2018.data.Sheet3(:,4)-FishRI18/1000).^2+(RI2019.data.Sheet3(:,4)-FishRI19/1000).^2+(RI2020.data.Sheet3(:,4)-FishRI20/1000).^2)/11;
MSPEFishC_RI=((RI2010.data.Sheet4(:,4)-FishRI10/1000).^2+(RI2011.data.Sheet4(:,4)-FishRI11/1000).^2+(RI2012.data.Sheet4(:,4)-FishRI12/1000).^2+(RI2013.data.Sheet4(:,4)-FishRI13/1000).^2+(RI2014.data.Sheet4(:,4)-FishRI14/1000).^2+(RI2015.data.Sheet4(:,4)-FishRI15/1000).^2+(RI2016.data.Sheet4(:,4)-FishRI16/1000).^2+(RI2017.data.Sheet4(:,4)-FishRI17/1000).^2+(RI2018.data.Sheet4(:,4)-FishRI18/1000).^2+(RI2019.data.Sheet4(:,4)-FishRI19/1000).^2+(RI2020.data.Sheet4(:,4)-FishRI20/1000).^2)/11;

MSPEFishA_VA=((VA2010.data.Sheet1(:,4)-FishVA10/1000).^2+(VA2011.data.Sheet1(:,4)-FishVA11/1000).^2+(VA2012.data.Sheet1(:,4)-FishVA12/1000).^2+(VA2013.data.Sheet1(:,4)-FishVA13/1000).^2+(VA2014.data.Sheet1(:,4)-FishVA14/1000).^2+(VA2015.data.Sheet1(:,4)-FishVA15/1000).^2+(VA2016.data.Sheet1(:,4)-FishVA16/1000).^2+(VA2017.data.Sheet1(:,4)-FishVA17/1000).^2+(VA2018.data.Sheet1(:,4)-FishVA18/1000).^2+(VA2019.data.Sheet1(:,4)-FishVA19/1000).^2+(VA2020.data.Sheet1(:,4)-FishVA20/1000).^2)/11;
MSPEFishM_VA=((VA2010.data.Sheet2(:,4)-FishVA10/1000).^2+(VA2011.data.Sheet2(:,4)-FishVA11/1000).^2+(VA2012.data.Sheet2(:,4)-FishVA12/1000).^2+(VA2013.data.Sheet2(:,4)-FishVA13/1000).^2+(VA2014.data.Sheet2(:,4)-FishVA14/1000).^2+(VA2015.data.Sheet2(:,4)-FishVA15/1000).^2+(VA2016.data.Sheet2(:,4)-FishVA16/1000).^2+(VA2017.data.Sheet2(:,4)-FishVA17/1000).^2+(VA2018.data.Sheet2(:,4)-FishVA18/1000).^2+(VA2019.data.Sheet2(:,4)-FishVA19/1000).^2+(VA2020.data.Sheet2(:,4)-FishVA20/1000).^2)/11;
MSPEFishL_VA=((VA2010.data.Sheet3(:,4)-FishVA10/1000).^2+(VA2011.data.Sheet3(:,4)-FishVA11/1000).^2+(VA2012.data.Sheet3(:,4)-FishVA12/1000).^2+(VA2013.data.Sheet3(:,4)-FishVA13/1000).^2+(VA2014.data.Sheet3(:,4)-FishVA14/1000).^2+(VA2015.data.Sheet3(:,4)-FishVA15/1000).^2+(VA2016.data.Sheet3(:,4)-FishVA16/1000).^2+(VA2017.data.Sheet3(:,4)-FishVA17/1000).^2+(VA2018.data.Sheet3(:,4)-FishVA18/1000).^2+(VA2019.data.Sheet3(:,4)-FishVA19/1000).^2+(VA2020.data.Sheet3(:,4)-FishVA20/1000).^2)/11;
MSPEFishC_VA=((VA2010.data.Sheet4(:,4)-FishVA10/1000).^2+(VA2011.data.Sheet4(:,4)-FishVA11/1000).^2+(VA2012.data.Sheet4(:,4)-FishVA12/1000).^2+(VA2013.data.Sheet4(:,4)-FishVA13/1000).^2+(VA2014.data.Sheet4(:,4)-FishVA14/1000).^2+(VA2015.data.Sheet4(:,4)-FishVA15/1000).^2+(VA2016.data.Sheet4(:,4)-FishVA16/1000).^2+(VA2017.data.Sheet4(:,4)-FishVA17/1000).^2+(VA2018.data.Sheet4(:,4)-FishVA18/1000).^2+(VA2019.data.Sheet4(:,4)-FishVA19/1000).^2+(VA2020.data.Sheet4(:,4)-FishVA20/1000).^2)/11;

subtable4Trips=vertcat(horzcat("Trips","alpha-maxmin","mlogit","clogit","CARA"),horzcat(vertcat("MA","RI","NY","CT","NJ","DE","MD","VA","NC"),100*vertcat(horzcat(mean(MSPEA_MA),mean(MSPEM_MA),mean(MSPEL_MA),mean(MSPEC_MA))/min(horzcat(mean(MSPEA_MA),mean(MSPEM_MA),mean(MSPEL_MA),mean(MSPEC_MA))),horzcat(mean(MSPEA_RI),mean(MSPEM_RI),mean(MSPEL_RI),mean(MSPEC_RI))/min(horzcat(mean(MSPEA_RI),mean(MSPEM_RI),mean(MSPEL_RI),mean(MSPEC_RI))),horzcat(mean(MSPEA_NY),mean(MSPEM_NY),mean(MSPEL_NY),mean(MSPEC_NY))/min(horzcat(mean(MSPEA_NY),mean(MSPEM_NY),mean(MSPEL_NY),mean(MSPEC_NY))),horzcat(mean(MSPEA_CT),mean(MSPEM_CT),mean(MSPEL_CT),mean(MSPEC_CT))/min(horzcat(mean(MSPEA_CT),mean(MSPEM_CT),mean(MSPEL_CT),mean(MSPEC_CT))),horzcat(mean(MSPEA_NJ),mean(MSPEM_NJ),mean(MSPEL_NJ),mean(MSPEC_NJ))/min(horzcat(mean(MSPEA_NJ),mean(MSPEM_NJ),mean(MSPEL_NJ),mean(MSPEC_NJ))),horzcat(mean(MSPEA_DE),mean(MSPEM_DE),mean(MSPEL_DE),mean(MSPEC_DE))/min(horzcat(mean(MSPEA_DE),mean(MSPEM_DE),mean(MSPEL_DE),mean(MSPEC_DE))),horzcat(mean(MSPEA_MD),mean(MSPEM_MD),mean(MSPEL_MD),mean(MSPEC_MD))/min(horzcat(mean(MSPEA_MD),mean(MSPEM_MD),mean(MSPEL_MD),mean(MSPEC_MD))),horzcat(mean(MSPEA_VA),mean(MSPEM_VA),mean(MSPEL_VA),mean(MSPEC_VA))/min(horzcat(mean(MSPEA_VA),mean(MSPEM_VA),mean(MSPEL_VA),mean(MSPEC_VA))),horzcat(mean(MSPEA_NC),mean(MSPEM_NC),mean(MSPEL_NC),mean(MSPEC_NC))/min(horzcat(mean(MSPEA_NC),mean(MSPEM_NC),mean(MSPEL_NC),mean(MSPEC_NC))))));
subtable4Fish=vertcat(horzcat("Landings","alpha-maxmin","mlogit","clogit","CARA"),horzcat(vertcat("MA","RI","NY","CT","NJ","DE","MD","VA","NC"),100*vertcat(horzcat(mean(MSPEFishA_MA),mean(MSPEFishM_MA),mean(MSPEFishL_MA),mean(MSPEFishC_MA))/min(horzcat(mean(MSPEFishA_MA),mean(MSPEFishM_MA),mean(MSPEFishL_MA),mean(MSPEFishC_MA))),horzcat(mean(MSPEFishA_RI),mean(MSPEFishM_RI),mean(MSPEFishL_RI),mean(MSPEFishC_RI))/min(horzcat(mean(MSPEFishA_RI),mean(MSPEFishM_RI),mean(MSPEFishL_RI),mean(MSPEFishC_RI))),horzcat(mean(MSPEFishA_NY),mean(MSPEFishM_NY),mean(MSPEFishL_NY),mean(MSPEFishC_NY))/min(horzcat(mean(MSPEFishA_NY),mean(MSPEFishM_NY),mean(MSPEFishL_NY),mean(MSPEFishC_NY))),horzcat(mean(MSPEFishA_CT),mean(MSPEFishM_CT),mean(MSPEFishL_CT),mean(MSPEFishC_CT))/min(horzcat(mean(MSPEFishA_CT),mean(MSPEFishM_CT),mean(MSPEFishL_CT),mean(MSPEFishC_CT))),horzcat(mean(MSPEFishA_NJ),mean(MSPEFishM_NJ),mean(MSPEFishL_NJ),mean(MSPEFishC_NJ))/min(horzcat(mean(MSPEFishA_NJ),mean(MSPEFishM_NJ),mean(MSPEFishL_NJ),mean(MSPEFishC_NJ))),horzcat(mean(MSPEFishA_DE),mean(MSPEFishM_DE),mean(MSPEFishL_DE),mean(MSPEFishC_DE))/min(horzcat(mean(MSPEFishA_DE),mean(MSPEFishM_DE),mean(MSPEFishL_DE),mean(MSPEFishC_DE))),horzcat(mean(MSPEFishA_MD),mean(MSPEFishM_MD),mean(MSPEFishL_MD),mean(MSPEFishC_MD))/min(horzcat(mean(MSPEFishA_MD),mean(MSPEFishM_MD),mean(MSPEFishL_MD),mean(MSPEFishC_MD))),horzcat(mean(MSPEFishA_VA),mean(MSPEFishM_VA),mean(MSPEFishL_VA),mean(MSPEFishC_VA))/min(horzcat(mean(MSPEFishA_VA),mean(MSPEFishM_VA),mean(MSPEFishL_VA),mean(MSPEFishC_VA))),horzcat(mean(MSPEFishA_NC),mean(MSPEFishM_NC),mean(MSPEFishL_NC),mean(MSPEFishC_NC))/min(horzcat(mean(MSPEFishA_NC),mean(MSPEFishM_NC),mean(MSPEFishL_NC),mean(MSPEFishC_NC))))));
Table4=vertcat(subtable4Trips,subtable4Fish);
clearvars subtable4Trips subtable4Fish

toc;

